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EXTRAGALACTIC MILLIMETER- WAVE SOURCES IN SOUTH POLE TELESCOPE SURVEY DATA: 
SOURCE COUNTS, CATALOG, AND STATISTICS FOR AN 87 SQUARE-DEGREE FIELD 
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ABSTRACT 

We report the results of an 87 deg 2 point-source survey centered at R.A. 5 h 30 m , decl. —55° taken 
with the South Pole Telescope (SPT) at 1.4 and 2.0 mm wavelengths with arc-minute resolution and 
milli-Jansky depth. Based on the ratio of flux in the two bands, we separate the detected sources 
into two populations, one consistent with synchrotron emission from active galactic nuclei (AGN) and 
one consistent with thermal emission from dust. We present source counts for each population from 
11 to 640 mJy at 1.4 mm and from 4.4 to 800 mJy at 2.0 mm. The 2.0 mm counts are dominated 
by synchrotron-dominated sources across our reported flux range; the 1.4 mm counts are dominated 
by synchroton-dominated sources above ~ 15 mJy and by dust-dominated sources below that flux 
level. We detect 141 synchrotron-dominated sources and 47 dust-dominated sources at S/N > 4.5 
in at least one band. All of the most significantly detected members of the synchrotron-dominated 
population are associated with sources in previously published radio catalogs. Some of the dust- 
dominated sources arc associated with nearby (z <C 1) galaxies whose dust emission is also detected 
by the Infrared Astronomy Satellite (IRAS). However, most of the bright, dust-dominated sources 
have no counterparts in any existing catalogs. We argue that these sources represent the rarest and 
brightest members of the population commonly referred to as submillimetcr galaxies (SMGs). Because 
these sources are selected at longer wavelengths than in typical SMG surveys, they are expected to 
have a higher mean redshift distribution and may provide a new window on galaxy formation in the 
early universe. 

Subject headings: galaxies: high-redshift — galaxies: surveys — submillimeter 
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1. INTRODUCTION 

The 10-meter South Pole Telescope (SPT, Carlstrom 
et al. 2009) is a millimeter/submillimetcr (mm/submm) 
telescope located at the geographic South Pole and de- 
signed for low-noise observations of diffuse, low-contrast 
sources such as anisotropy in the cosmic microwave back- 
ground (CMB). The first camera installed on the SPT is 
a 960-element bolomctric receiver designed to perform a 
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mass-limited survey of galaxy clusters via their Sunyaev- 
Zel'dovich (SZ) signature (Sunyaev & Zeldovich 1972) 
over a large area of the southern sky. This survey is cur- 
rently underway, and the SPT team recently published 
the first-ever discovery of galaxy clusters through their 
SZ signature (Staniszcwski ct al. 2009, hereafter S09). 

The sensitivity and angular resolution of the SPT make 
it an excellent instrument for detecting extragalactic 
sources of emission. In this work, we report on source 
detections in a small part of the SPT survey, namely a 
single 87 deg 2 field centered at right ascension (R.A.) 
5 h 30 m , declination (dccl.) -55° (J2000). This field was 
surveyed by the SPT in the 2008 season to roughly mJy 
depth at 1.4 mm and 2.0 mm (220 and 150 GHz). The 
data presented here represent a major step forward in 
mm source detection at these flux levels, both in area 
surveyed and in the ability to distinguish between source 
populations using internal estimates of source spectral 
properties. 

Simultaneous information in two bands for each de- 
tected source should allow us to separate our detections 
into distinct source populations. Based on previous sur- 
veys at mm wavelengths and on surveys in neighbor- 
ing centimeter (cm) and submm bands, we expect the 
sources we detect to fall into two broad categories: 1) 
sources with flat or decreasing brightness with decreasing 
wavelength, consistent with synchrotron emission from 
active galactic nuclei (AGN, typically S oc A"' 1 ); and 
2) sources with increasing brightness with decreasing 
wavelength, consistent with thermal emission (typically 
S oc A~ -3 ) from dust-enshrouded star-forming galaxies. 

The synchrotron-dominated source population is well- 
established from radio surveys (see De Zotti et al. (2010) 
for a recent review). Despite these sources' decreasing 
brightness from radio to mm wavelengths, simple extrap- 
olations of radio and cm counts of these sources to the 
SPT bands predict that we should detect a significant 
number of these sources. This prediction is bolstered by 
the results of 3 mm follow-up of 1.5 cm-detected sources 
presented in Sadler ct al. (2008), which showed that these 
sources still emitting strongly at mm wavelengths, and 
by detections of synchrotron-dominated AGN emission 
made in mm/submm surveys much smaller than the SPT 
survey but at similar depths (Voss et al. 2006). Measure- 
ments of the mm fluxes of a large sample of these sources 
have the power to inform both astrophysical models of 
their emission and predictions for the extent of their con- 
tamination to the CMB power spectrum (e.g., Toffolatti 
ct al. 2005; Reichardt et al. 2009) and the SZ signal from 
galaxy clusters (e.g., Lin & Mohr 2007; Sehgal et al. 
2010). 

The dust-enshrouded star-forming galaxy population 
has been the subject of considerable interest for over 
thirty years (sec Rickc & Lcbofsky (1979) for a pre- 
IRAS review), but the Infrared Astronomy Satellite 
(IRAS) was the first instrument to systematically dis- 
cover such objects (Sanders & Mirabel 1996). Optical 
and UV observations show these sources to be heavily 
dust-obscured, and to typically have disturbed morpholo- 
gies and high star formation rates, indicative of recent 
or ongoing mergers (Lagache et al. 2005). The emis- 
sion from IR to mm wavelengths in these sources arises 
from short-wavelength photons emitted by stars, which 



are absorbed by dust grains and re-radiated at longer 
wavelengths (Draine 2003). 

Measurements of the cosmic infrared background 
(CIB) show that over half the energy emitted since the 
surface of last scattering has been absorbed and re- 
radiated by dust (Dwck et al. 1998). IRAS, however, de- 
tected mostly low-rcdshift (z < 1) objects, and these rel- 
atively rare and nearby sources contributed only a small 
fraction of the CIB (Lc Floc'h et al. 2005; Caputi et al. 
2007), indicating that the bulk of CIB sources are at 
high redshift. The first systematic survey of high red- 
shift sources which contribute significantly to the CIB 
— the population now known as submillimeter galaxies 
(SMGs) — was carried out a decade ago at 850 fxm by 
the SCUBA instrument on the 15-m JCMT telescope 
(Smail et al. 1997; Hughes et al. 1998; Barger et al. 1998; 
Holland et al. 1999). Owing to the spectrum of SMGs — 
a modified ~ 30 K blackbody that rises steeply with 
decreasing wavelength, counteracting the expected flux 
diminution with redshift — they can be detected indepen- 
dently of redshift from roughly 500 fim to 2 mm (Blain 
et al. 2002). This implies that the source luminosity is 
roughly proportional to the brightness from 1 < z < 10. 

Hundreds of SMGs have now been detected by ground- 
based telescopes in surveys of blank fields, but only over 
a total area on the order of a square degree (Scott et al. 
2002; Borys et al. 2003; Greve et al. 2004; Laurent et al. 
2005; Coppin et al. 2006; Bertoldi et al. 2007; Pereraet al. 
2008; Scott et al. 2008; Austermann et al. 2010; Weifi 
et al. 2009). Recently, results were published from the 
Balloon-borne Large-Aperture Submillimeter Telescope 
(BLAST) which surveyed nearly ten square degrees at 
250, 350, and 550 /im and measured important proper- 
ties such as dust temperatures and clustering amplitude 
for SMGs (Devlin et al. 2009; Patanchon ct al. 2009; Dye 
et al. 2009; Viero et al. 2009). The discovery and study 
of SMGs has revolutionized our understanding of galaxy 
formation by providing a view of galaxy formation which 
is both unbiased with redshift and inaccessible to optical 
surveys. Observations of these objects (see Blain et al. 
(2002) for a review) indicate that: 1) they have dynam- 
ical masses of ~10 11 M Q and total far-infrared luminosi- 
ties of-10 13 L Q (Swinbank et al. 2004; Greve et al. 2005; 
Chapman et al. 2005; Kovacs et al. 2006; Pope et al. 
2006); 2) they are forming stars prodigiously at ~ 1000 
M /ycar (Chapman ct al. 2005; Tacconi et al. 2006); 3) 
their abundance appears to peak at z ~ 2.5 (Pope et al. 
2005; Chapman et al. 2005; Aretxaga et al. 2007; Chapin 
et al. 2009); 4) from observations (Tacconi et al. 2008) 
and simulation (Barnes & Hernquist 1991; Narayanan 
ct al. 2010), the prodigious star formation rate seen in 
SMGs is believed to be intrinsically linked to mergers; 
and 5) SMGs are an early phase in the formation of the 
most massive galaxies and are among the largest gravi- 
tationally collapsed objects in this early epoch of galaxy 
formation (Blain et al. 2004; Swinbank et al. 2008). 

In the context of this broad field of IR/submm/mm- 
selected, dust-enshrouded, star-forming galaxies, this pa- 
per presents the detection in the SPT data of a pop- 
ulation of dust-dominated sources with surprising and 
intriguing properties. These sources are significantly 
brighter and rarer than the submm-selected population 
in the literature. Furthermore, the majority of these 
sources do not have counterparts in IRAS, indicating 
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that they are not members of the standard local ultra- 
luminous infrared galaxy (ULIRG) population. This ap- 
parently new family of sources represents the most sig- 
nificant new result of this work. 

This paper is divided into several sections. Sec. 2 de- 
scribes the SPT observations, data reduction, matched 
filter for point-source signal, and source-finding algo- 
rithm. (The observations and data reduction through 
the mapmaking step are described in greater detail in 
S09 and Carlstrom et al. 2009.) Sec. 3 discusses the 
properties of the filtered maps, presents the source cat- 
alog, describes our procedures for checking astrometry 
and estimating completeness and purity, and discusses 
basic source properties, including raw spectral classifi- 
cation. Sec. 4 describes our procedure for estimating 
each source's intrinsic flux and spectral index, which 
we use to separate our sources into two spectrally dis- 
tinct (synchrotron-dominated and dust-dominated) pop- 
ulations. The statistical method used for flux estimation 
is described in detail in a companion paper, Crawford 
et al. (2010). Sec. 5 presents source counts for each band. 
Sec. 6 discusses associations with external catalogs. Fi- 
nally, Sec. 7 presents counts for each of the populations 
and discusses the implications, including the potential 
for a newly discovered population of sources. 

In a companion paper (Hall et al. 2010), we present the 
spatial power spectra of the sources below our detection 
threshold. 

2. OBSERVATIONS, DATA REDUCTION, AND SOURCE 
FINDING 

This work is based on observations of a single ~ 100- 
deg 2 field. The timestream data for each observation, 
constituting a single pass over the field, were processed 
and combined to make a map of the field for each ob- 
serving band. The maps from several hundred individual 
observations of the field were combined and converted to 
CMB fluctuation temperature units using a calibration 
from the CMB anisotropy as measured by the Wilkinson 
Microwave Anisotropy Probe (WMAP, Hinshaw et al. 
(2009)). Each single-band map is filtered to optimize 
point-source detection. A variant of the CLEAN algo- 
rithm (Hbgbom 1974) was then used to search for sources 
in an 87 deg 2 sub-region of the filtered maps. Finally, 
amplitudes of the detected peaks in the filtered maps 
were converted from CMB fluctuation temperature units 
to flux (in units of Jy). 

2.1. Observations 

During the 2008 observing season, the 960-element 
SPT camera included detectors sensitive to radiation 
within bands centered at approximately 1.4 mm, 2.0 mm, 
and 3.2 mm (220 GHz, 150 GHz, and 95 GHz). The first 
field mapped to the targeted survey depth was centered 
at R.A. 5 h 30 m , deck -55° (J2000). Results in this paper 
are based on 607 hours of observing time, using only the 
1.4 mm and 2.0 mm data from the 87 deg 2 portion of 
the field that was mapped with near-uniform coverage. 
(Every 1 arcmin patch in the included area is required to 
have uniform coverage to 10%.) The scanning strategy 
consisted of constant-elevation scans across the 10° wide 
field. After each single scan back and forth across the 
field, the telescope executes a 0.125° step in elevation. 
A complete set of scans covering the entire field takes 



approximately two hours, and we refer to each complete 
set as an observation. Between individual observations 
of the field, we perform a series of short calibration mea- 
surements, including measurements of a chopped ther- 
mal source, 2 degree elevation nods, and scans across the 
galactic HII regions RCW38 and MAT5a. This series 
of regular calibration measurements allows us to identify 
detectors with good performance, assess relative detector 
gains, monitor atmospheric opacity and beam parame- 
ters, and constrain pointing variations. S09 and Carl- 
strom et al. (2009) describe details of the observations of 
this field and the telescope. 

2.2. Data Reduction 

The reduction of SPT data for this work up to and 
including the mapmaking step is very similar to the re- 
duction pipeline used to produce the maps used in S09, 
where details of the analysis can be found. Broadly, the 
pipeline consists of filtering the time-ordered data from 
each individual detector, reconstructing the pointing for 
each detector, and combining data from all detectors in 
a given observing band into a map by simple inverse- 
variance-weighted binning and averaging. 

The filtering used in this analysis differ from that in 
S09. In this work the time-ordered detector data were 
filtered with a 0.18 Hz Fourier-domain high-pass filter. 
With our scan speeds, the high pass filter removes spa- 
tial scales > 45'. We project out a common mode which 
consists of three spatial modes (mean, and tilts along 
two axes) constructed from the mean of all working de- 
tectors in a single band, weighted by the x and y position 
in the focal plane. As atmospheric signal is highly cor- 
related between detectors, removing this common mode 
should eliminate the majority of the atmospheric fluctu- 
ation power in the detector timestreams. The common- 
mode subtraction acts as a spatial high-pass filter with a 
characteristic scale that roughly corresponds to the one 
degree angular field of view of the array. This filter op- 
tion was demonstrated to remove more atmosphere from 
the timestream than the method described in S09, but 
its choice was not critical. As the common mode is con- 
structed independently for each band, the response to 
large spatial modes on the sky can be slightly different 
between bands, but these modes are heavily de-weighted 
in the filter we later apply to the maps to enhance point- 
source signal-to-noise. Finally, in contrast to S09, we did 
not mask bright point sources when filtering the time- 
ordered data, because we wanted to ensure that the filter 
transfer function would be the same for all sources in the 
maps. 

2.3. Flux Calibration and Beam Measurements 

The relative gains of the detectors and their gain vari- 
ations over time were estimated using measurements of 
their response to a chopped thermal source. (These mea- 
surements took place before each observation of the sur- 
vey field, or every 2 hours.) This was the same relative 
calibration method used in S09. The absolute calibra- 
tion is based on the comparison of SPT measurements of 
degree-scale CMB fluctuations at 2.0 mm directly to the 
WMAP 5-year maps. This was done using short, dedi- 
cated observations of four large fields, totaling 1250 deg 2 . 
Details of the cross-calibration with WMAP are given in 
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Lueker et al. (2010). We estimate the uncertainty of this 
calibration to be 3.6%. We applied this calibration to 
our 1.4 mm band by comparing 2.0 mm and 1.4 mm es- 
timates of CMB anisotropy in our deep survey regions. 
This internal cross-calibration for 1.4 mm is consistent 
with a direct absolute calibration from RCW38, but with 
higher precision. We estimate the 1.4 mm calibration 
uncertainty to be 7.2%. Because the 1.4 mm calibration 
is derived from the 2.0 mm calibration to WMAP, the 
calibration uncertainties in the two bands are correlated 
with an estimated correlation coefficient of roughly 0.5. 

Main-lobe beams were measured using the brightest 
sources in the field and are adequately fit by 2d Gaus- 
sians with full width at half-maximum (FWHM) equal 
to 1.05' and 1.15' at 1.4 mm and 2.0 mm. Large-angle 
sidelobes were measured using planet observations, but 
the angular scales on which these sidelobes are important 
are heavily downweightcd in the filter. We estimate that 
beam-shape uncertainties contribute roughly 2% and 5% 
to the absolute flux estimates in our 2.0 mm and 1.4 mm 
bands. This uncertainty was added in quadrature to the 
calibration uncertainty in our flux estimates. A visual 
representation of the SPT 2.0 mm beam is shown in 
Padin et al. (2008), and symmetrized ^-space profiles for 
the beams at both wavelengths are shown in Lueker et al. 
(2010). 

A subtlety in estimating the spectral index (Eq. 5, be- 
low) is that the effective band centers (which fold into the 
index determination) depend upon spectral index. Using 
the measured passbands for 1.4 mm and 2.0 mm, we find 
that if one were to assume an index a = — 1 in the deter- 
mination of the band centers, a source with a — 3 would 
be measured with a 2% bias in the spectral index. In ad- 
dition, the beam shape (and so flux) will change with the 
spectral index. These subtleties can both be neglected 
to the accuracy of the results presented here. 

2.4. Source Extraction 

2.4.1. Matched Filter 

We enhanced the point-source signal-to-noise ratio in 
the SPT maps by applying a matched spatial filter (see 
e.g., Tegmark & de Oliveira-Costa 1998) to each single- 
band map. The matched filter combines knowledge of 
the instrument beam and any other filtering that has 
been performed on the data with an estimate of noise 
covariance to optimize the signal-to-noise of a source in 
the filtered map. This matched filter ip is applied in the 
Fourier domain and is given by: 

i, = —== 1 

V^N- 1 T 

where N is the noise covariance matrix (includ- 
ing astrophysical contaminants such as primary CMB 
anisotropy), and t is the assumed source shape in the 
map, which in the case of point sources is a function of 
beam and filtering only. Real-space visual representa- 
tions of r for 1.4 mm and 2.0 mm — with filtering very 
similar to that used in the maps in this work — are shown 
in Plagge et al. (2010). 

As in S09, the instrumental and atmospheric contribu- 
tions to the noise covariance in each band were estimated 
by computing the average power spectrum of hundreds 
of signal-free maps, constructed by multiplying half of 



the 320 individual-observation maps of the field by -1, 
half by +1, and then summing. (Hereafter, this will be 
referred to as the "difference" map, in which all astro- 
physical signal has been removed from the map, but the 
atmospheric and detector noise remains.) The main as- 
trophysical contribution to the noise covariance is ex- 
pected to be primary CMB anisotropy, so an estimate of 
the CMB power spectrum was added to the noise covari- 
ance. Adding further astrophysical contributions such as 
the SZ background and point sources below our detection 
threshold has a negligible effect on our results. 

The source shape used in the matched filter is the 
convolution of our measured beam and the map-domain 
equivalent of any timestream filtering we have performed. 
We measured the effect of timestream filtering on the ex- 
pected shape of point sources in our maps by perform- 
ing signal-only simulations of our data processing, using 
a single model point source projected into our detector 
timestreams as the only timestream component. 

Observations used here were performed using constant 
declination scans. The maps were pixelized using a flat- 
sky projection of the sphere in which the mapping of 
right ascension to map rows is a function of position in 
the map. 27 As a result, the effects of timestream filtering 
on source shape are also map-position-dependent. To ac- 
count for this, we break the single-band coadded signal 
maps into nine tiles and perform our signal-only simula- 
tions nine times — once with the model source located at 
the center of each tile. We broke the map into nine tiles 
(as opposed to four or sixteen) as it solved the problem 
with the greatest economy. We estimated the noise co- 
variance separately for each tile, as the projection of non- 
white timestream noise into the map is also a function of 
position. The noise in the outer tiles is within 5% of the 
noise in the center tile (see Sec. 3). We constructed nine 
matched filters from these inputs and performed source 
finding on each map tile individually with the matched 
filter constructed from that tile's inputs. 

2.4.2. Candidate Identification 

Source candidates were identified in the filtered maps 
using a variant of the CLEAN algorithm from radio as- 
tronomy (Hogbom 1974). Briefly, the CLEAN proce- 
dure involves iteratively identifying the highest peak in 
the filtered map and subtracting off a model for the 
source shape centered on that peak until no peaks are 
left above the detection threshold. To account for several 
non-idealities, including finite-sized map pixels, slightly 
imperfect source shape models, and possibly extended 
sources, the source model subtraction is performed with 
a multiplicative factor less than unity, usually called the 
loop gain (after the analogous parameter in electronic 
feedback circuits). We use a loop gain of 0.1 in this 
work. 

In interferometric radio observations, the source shape 
template in the CLEAN process is the interferometer's 
"dirty beam" ; in our case, it is the source shape in the 
filtered maps: 

t' = V> t. (2) 
As discussed in the previous section, the matched filter 

27 We chose this pixelization because it minimizes beam distor- 
tions, which are significant in flat-sky pixclizations in which pixel 
rows arc at constant declination. 
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ip was independently calculated for nine different region 
of each band's map in order to account for the map- 
position-dependent shape of the noise and filtering. In 
constructing the source shape template ipr, we used the 
appropriate version of -0 depending on the position of 
the peak being CLEANcd. For each source peak we in- 
dependently calculated a pre-matched-filter source shape 
r to account for this positional dependence. In the map 
pixelization used, this calculation consists of a simple ro- 
tation of a fiducial source shape, so this step was not 
unduly computationally intensive. 

We ran our version of CLEAN on each band's filtered 
map individually until there were no peaks above 3a 
left in the map. All map pixels identified above the 3cr 
threshold were then sorted by significance and gathered 
into discrete sources using an association radius between 
30 arcsec and 2 arcmin, depending on the brightness of 
the source. In other words, the brightest pixel found by 
CLEAN was declared to be the first source, then the list 
of pixels was examined in descending order of brightness, 
asking if each pixel should be declared a new source or 
associated with a source already identified. Source fluxes 
were assigned by converting the value in the filtered map 
of the brightest pixel associated with a source from CMB 
fluctuation temperature units to flux (in units of Jy) us- 
ing the following relation: 

ott i rn a <~i i n26 ( ^b^cmb \ x^e x 
S[Jy] = T pcak x Afi f x 10 26 x — — J (e x_ 1)2 - 

(3) 

where x = ^/(^bJcmb) and Afl{ is defined by: 



dudv ip(u,v) t(u,v) 



(4) 



which can be thought of as the effective solid angle under 
the filtered source template, in that a point source of flux 
S will have peak brightness S/AQf in the filtered map. 
This flux estimate will be biased for extended sources, 
which are discussed in Sec. 3.5. Source positions were 
obtained by calculating the center of brightness of all 
pixels (each pixel being 0.25' x 0.25') associated with a 
given source. 

3. MAPS AND CATALOG 

The filtered 2.0 mm and 1.4 mm maps used for source 
candidate identification are shown in Fig. 1 and Fig. 2. 
The total area shown in each map is 87 square degrees. 
As previously noted, the noise varies at the level of ±5% 
across the maps, mainly as a function of declination. 
This trend with declination is due to the fact that the 
coverage is nearly uniform in right ascension, resulting in 
coverage per unit solid angle that varies as cos(dccl.) (i.e. 
the noise is systematically 5% lower at decl.=— 60 than 
at dccl.=— 55). The typical RMS of the map is 1.3 mJy 
at 2.0 mm and 3.4 mJy at 1.4 mm. The noise distribu- 
tion closely approximates a Gaussian, as is evident from 
the central part of the pixel distributions shown in Fig. 3. 
The fact that the maps are so uniform and the noise is 
so well-understood makes the analysis much easier and 
gives us confidence in the robustness of our results. 

Detections in both bands are listed in the final catalog 
as a single source if they are offset < 30 arcsec between 
bands. For sources detected in both bands, we adopt the 



position of the more significant detection. We are far 
enough above the confusion limit that this simple and 
intuitive method is adequate. For sources detected in 
only one band we use the flux in the cleaned map for the 
second band at the position of the detection. Table 5 lists 
the properties of every source candidate detected above 
4.5(7 in either the 2.0 mm or 1.4 mm band, while the full 
table listing all 3496 sources above 3(7 in cither map will 
be available in machine-readable form in the electronic 
version of this article; for now, the full list is available 
from the public SPT website. 28 

3.1. Catalog field descriptions 

Each column in Table 5 corresponds to a field in the 
87 deg 2 SPT two-band source catalog. Descriptions of 
the catalog fields / table columns are as follows: 

1. Source ID: the IAU designation for the SPT- 
detected source. 

2. RA: right ascension (J2000) in degrees. 

3. DEC: declination (J2000) in degrees. 

4. S/N (2.0 mm): detection significance (signal-to- 
noise ratio) in the 2.0 mm band. 

5. <; raw (2.0 mm): raw flux (uncorrected for flux 
boosting) in the 2.0 mm band. 

6. S dlst (2.0 mm): de-boosted flux values encompass- 
ing 16%, 50%, and 84% (68% probability enclosed, 
or lcr for the equivalent normal distribution) of 
the cumulative posterior probability distribution 
for 2.0 mm flux, as estimated using the procedure 
described in Sec. 4. 

7. S/N (1.4 mm): same as (4), but for 1.4 mm. 

8. S"' aw (1.4 mm): same as (5), but for 1.4 mm. 

9. 5 dlst (1.4 mm): same as (6), but for 1.4 mm. 

10. a raw : estimate (from the raw flux in each band) of 
the 2.0 mm— 1.4 mm spectral index a, where a is 
the slope of the (assumed) power-law behavior of 
source flux function of wavelength: 



S oc X a . 



(5) 



11. a dist : 16%, 50%, and 84% estimates of the spec- 
tral index, based on the probability distributions 
for spectral index estimated using the procedure 
described in Sec. 4. 

12. P(a > 1.66): fraction of the spectral index pos- 
terior probability distribution above the threshold 
value of 1.66. A higher value of P means the source 
is more likely to be dust-dominated. 

13. Type: source classification (synchrotron- or dust- 
dominated), based on whether P(a > 1.66) is 
greater than or less than 0.5. 

28 http:/ /pole. uchicago.edu/public/data/vieira09/ 
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Fig. 1. — The filtered 2.0 mm map in a flat sky projection. The total sky area is 87 deg 2 and the field center is right ascension 
R.A. 5 h 30 m , deck -55° (J2000). The map is oriented such that north (increasing deck) is up and east (increasing R.A.) is 
left. Each pixel is 0.25' x 0.25'. The RMS in the map is 1.3 mjy, and the gray scale is ±5<r; the brightest source is > 500<r 
(> 650 mjy), and the scale saturates for most of the sources visible here. Because of time domain filtering, the source signal 
produces an arc from the impulse response of the filter as the detectors scan left and right across the field. The azimuthally 
symmetric ringing around bright sources is due to spatial high-pass filtering in the point-source matched filter. 
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14. Nearest SUMSS source: angular distance (in ar- 
seconds) from the nearest source in the 36 cm 
(843 MHz) Sydney University Molongolo Sky Sur- 
vey (SUMSS) (Mauch et al. 2003). There are 2731 
SUMSS sources in the SPT survey area. For a 1 ar- 
cmin association radius there is a 2.7% chance of 
random association for each SPT source. 

15. Nearest RASS source: angular distance (in arsec- 
onds) from the nearest source in the ROSAT All- 
Sky Survey (RASS) Bright Source Catalog (Voges 
ct al. 1999) or Faint Source Catalog (Voges et al. 
2000). There arc 1441 RASS sources in the SPT 
survey area. For a 1 arcmin association radius there 
is a 1.4% chance of random association for each 
SPT source. 

16. Nearest IRAS source: angular distance (in arsec- 
onds) from the nearest source in the IRAS Faint- 
Source Catalog (IRAS-FSC, Moshir ct al. 1992). 
There are 493 IRAS sources in the SPT survey 
area. For a 1 arcmin association radius there is 
a 0.8% chance of random association for each SPT 
source. 



3.2. Astrometry 

SPT pointing is reconstructed through a combination 
of an online pointing model (tied to regular observations 
with optical star cameras), corrections based on observa- 
tions of galactic HII regions (performed many times each 
observing day), and information from thermal and lin- 
ear displacement sensors on the telescope. The pointing 
reconstruction process is described in more detail in S09 
and Carlstrom et al. (2009). In this work we calibrate the 
absolute positions in the maps by comparing our best-fit 
positions for bright sources (S/N > 10) in our catalog 
with external determinations of those positions from the 
1.5 cm Australia Telescope 20 GHz Survey (AT20G) cat- 
alog, in which the absolute astrometry is tied to VLBI 
calibrators and is accurate at the 1-arcsec level (Murphy 
et al. 2010). We used 15 point sources for the 1.4 mm 
absolute astrometry correction, and 26 point sources for 
2.0 mm. Fig. 4 shows the distribution of offsets between 
SPT-dctcrmined positions and AT20G positions. 

3.3. Completeness and Purity 
3.3.1. Completeness 

We follow Scott et al. (2008) and estimate our com- 
pleteness by placing simulated sources in the actual sig- 
nal maps and performing the source extraction as with 
the real data. For the simulated source profile, we use 
the measured beam convolved with the map-domain es- 
timate of our timestream filtering and the matched filter. 
As with the matched filter and the CLEAN process, we 
use a different simulated source profile in each of the 
nine map tiles (see Sec. 2.4 for details). The simulated 
source is considered detected if it would have made it 
into our catalog — i.e., if it is detected by the source ex- 
traction algorithm at > 3er. As expected for maps whose 
variance is nearly uniform and is dominated by random, 
Gaussian-distributed noise, our cumulative completeness 
curves (fraction of simulated sources detected above a 



given flux) are fit well by error functions, as shown in 
Fig. 5. The exact functional form used here is 



/couiplv*^) 



1 



e-( s '~ s °) 2 / 2 ° 2 dS'. (6) 



On the basis of this test and the error-function fits, we 
expect the full > 3c catalog to be 50% complete at 3.5 
and 9.1 mJy in the in the 2.0 and 1.4 mm bands and to 
be 95% complete at 5.5 mJy and 14.1 mJy in the 2.0 and 
1.4 mm bands. 

3.3.2. Purity 

There is some ambiguity in the definition of "purity" 
or "false detection" when one is dealing with a source 
population with steep differential number counts, espe- 
cially if the detected fluxes are anywhere near the con- 
fusion limit. In such a situation, there will be at least 
one source at a non-negligible fraction of the detection 
threshold in every beam. In this work, we have chosen to 
define a false detection as a fluctuation above the detec- 
tion threshold in the absence of any mean point source 
contribution to the maps. We treat the problem of low- 
flux sources scattering above the detection threshold in 
the context of flux boosting in Sec. 4. 

We estimate our purity using two different methods, 
both of which are fairly common in the SMG literature 
(e.g., Perera et al. 2008). First, we invert our maps and 
run the matched filter and source-finding algorithm on 
the negative maps. This method is complicated by the 
fact that, at 2.0 mm, we expect to have real negative 
signal near the beam scale due to the thermal SZ signal 
from galaxy clusters. To deal with this, we mask the 
inverted 2.0 mm map around SZ cluster candidates de- 
tected at > 4.5cr. These candidates are identified using 
a filter optimized for extended sources with a particu- 
lar spatial profile (in this case a spherical (3 model, see 
S09 for details), so we should not be masking point-like 
noise fluctuations with this procedure. Our second esti- 
mate of purity comes from running the matched filter and 
source-finding algorithm on simulated maps. These sim- 
ulated maps contain atmospheric and instrumental noise 
(taken from our difference map - see Sec. 2.4), a real- 
ization of the CMB, and a white, Gaussian noise term 
meant to approximate the contribution from the back- 
ground of sources below the detection threshold. The 
results from both tests are shown in Fig. 6. In all cases, 
the quantity plotted is 



/p 



1 



N, 



False 



total 



(7) 



where iVf a i se is the number of false detections (as esti- 
mated, alternately, by one of the two methods described 
above) above a given S/N, and iV-totai is the total num- 
ber of detections above a given S/N in the real map. 
Both methods agree that at S/N > 4.5 our sample is 
> 90% pure. Perera ct al. (2008) argue that both of 
these methods will overestimate the true false detection 
rate, and this hypothesis is supported by the fraction of 
our synchrotron-dominated sources that have clear coun- 
terparts in other catalogs and/or our ATCA follow-up 
observations (see Sec. 6 for details). 

3.3.3. Contamination at 2.0 mm from SZ 
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Fig. 2. — The filtered 1.4 mm map. The RMS in the map is 3.4 mjy, and the gray scale is ±5<r; the brightest source is > 150<r 
(> 550 mjy). See the corresponding Fig. 1 for comments common to both maps. 
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Fig. 3. — The distribution of fluxes in map pixels. For each 
band, the lines are as follows: solid: The coadded signal map; 
dashed: The coadded difference map (see Sec. 2.4); dotted: 
Fit to the signal map pixel histogram noise peak. For each 
band, the fit is done to the full signal map and gives a = 
1.3 mJy for 2.0 mm and a = 3.4 mJy for 1.4 mm. The noise 
across the map is Gaussian. The negative tails are mainly 
due to ringing from the various effective high-pass filters on 
the sources in the map. 

In addition to complicating the purity analysis in the 
previous section, SZ decrements from galaxy clusters 
have the potential to contaminate our source measure- 
ments at 2.0 mm (though not at 1.4 mm, which is very 
close to the thermal SZ null). We believe that this con- 
tamination will be negligible at the source flux levels 
considered here for two reasons. One reason is because 
clusters are expected to be at least partially resolved by 
the SPT at 2.0 mm, meaning that their contribution to 
maps filtered to optimize point-source sensitivity will be 
diminished. The other reason is that the number den- 
sity of clusters with decrements deep enough to signifi- 
cantly affect the source fluxes presented here is expected 
to be quite low. The SZ contamination will be some- 
what boosted by the fact that the sources we investigate 
here are expected to be spatially correlated with galaxy 
clusters at some level (e.g., Coble et al. 2007; Bai et al. 



Fig. 4. — A comparison of relative pointing offsets between 
the SPT bands and AT20G catalog sources. Only sources 
with S/N > 10 which have a robust counterpart within 20" 
have been plotted. The errors for the SPT positions were 
estimated following Ivison et al. (2007). The errors plotted 
here are a quadrature sum of the SPT error and the quoted 
error from the AT20G catalog. The RMS for both SPT bands 
is < 3". 

2007), but we expect the net effect to be negligible even 
after accounting for this correlation. 

To make these arguments more quantitative, we inves- 
tigate the level of flux decrements in simulated SZ maps 
filtered in the same way the SPT data is filtered in this 
work. We take simulated 2.0 mm maps created using 
the technique described in Shaw et al. (2009), filter them 
with the estimate of the SPT 2.0 mm beam and filtering 
discussed in Sec. 2.4.1 to simulate SPT observation and 
data processing, and further filter them with the matched 
filter from Sec. 2.4.1. Finally, we convert the filtered map 
from temperature to flux (as in Sec. 2.4.2) and record the 
decrement in Jy at each simulated cluster location. We 
find roughly five clusters per square degree with at least 
a 1.3 mJy decrement in the filtered map — equivalent to 
a la noise fluctuation in the 2.0 mm SPT map. We find 
roughly one cluster per ten square degrees with at least 
a 5.8 mJy decrement in the filtered map — equivalent to 
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Fig. 5. — The left panel shows the results of the completeness simulation at 2.0 mm; the right panel shows the results of the 
completeness simulation at 1.4 mm. In each plot, the symbols with error bars show the fraction of sources recovered at > 3<r 
with error bars estimated from binomial statistics. The dashed line shows the best-fit model of the form shown in Eq. 6. 




Fig. 6. — Purity in the 2.0 mm-selected sample (left) and the 1.4 mm-selected sample (right). In each plot, the solid line 
indicates the purity (see Eq. 7) calculated using the inverted map to estimate the number of false detections, while the dashed 
line indicates the purity calculated using simulated maps to estimate the number of false detections. The situation is more 
complicated at 2.0 mm than at 1.4 mm due to the presence of SZ (see Sec. 3.3.2 for details). 
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a 4.5cr noise fluctuation in the 2.0 mm SPT map. 

The SPT beam is roughly 1 arcmin wide at 2.0 mm, 
meaning that there are over 1000 independent resolu- 
tion elements in one square degree of the SPT 2.0 mm 
map. Thus, if galaxy clusters were randomly distributed 
with respect to sources in the SPT maps, we would ex- 
pect fewer than 0.5% of sources to suffer ak or worse 
systematic flux reduction due to cluster SZ signal, and 
we would expect fewer than 0.01% of 4.5a sources to be 
completely canceled by an SZ decrement. Coble et al. 
(2007) estimate that radio sources are a factor of 9 ± 4 
times more likely to be found along a line of sight within 
0.5 arcmin of a cluster than in the field. This boosts 
the chance of a systematic la flux error to < 5% and 
the chance of a 4.5er error to < 0.1%. Dust-dominated 
sources are expected to be less correlated with galaxy 
clusters than radio sources Bai et al. (2007), so the ef- 
fect will be even smaller for these sources. The effect of 
this level of systematic flux error on our final counts is 
completely subdominant to the statistical uncertainty. 

3.4. Raw Spectral Classification and Source Association 

Based on previous surveys of sources at other wave- 
lengths, we expect most SPT sources to belong to one 
of two populations: one dominated by synchrotron emis- 
sion, the members of which should have an emission spec- 
trum that is flat or falling with decreasing wavelength, 
and one dominated by thermal emission of reprocessed 
starlight by dust, the members of which should have an 
emission spectrum that increases with decreasing wave- 
length. Our results confirm this picture. Of course, any 
individual source may have components of each in its 
emission, and the local slope of the spectral energy dis- 
tribution (SED) will be further modulated by the redshift 
of the source. Though our actual source characterization 
is based on the integrated posterior probability density 
function (PDF) of the spectral index, estimated using the 
method described in Sec. 4.2, a plot of raw 1.4 mm flux 
vs. raw 2.0 mm flux, as in Fig. 7, gives the basic picture. 
Of the sources detected above 4.5er in both bands, the 
synchrotron-dominated sources occupy a locus of points 
close to the line a = — 1, where the spectral index a 
is defined in Eqn. 5. The dust-dominated sources de- 
tected in both bands occupy a clearly separated locus of 
points close to the a = 3 line. Also worth noting in this 
plot is that effectively all of the high-S/N synchrotron- 
dominated sources have counterparts in external cata- 
logs, while many of the high-S/N dust-dominated sources 
do not. This point is explored in greater detail in Sec. 6 
and Sec. 7.2. 

3.5. Extended Sources and Other Notes 

As is evident from Eq. 3 and Eq. 4, our flux estimates 
rest on the assumption that all sources have the same 
shape in our filtered maps. Since the assumed source 
shape is just that of our beam and filtering, this as- 
sumption will only be valid for point-like objects. This 
method will not provide accurate flux estimates for re- 
solved sources. For example, our method will underesti- 
mate the flux of a source with a FWHM = 0.25 arcmin 
Gaussian profile by 3% at 2.0 mm and 4% at 1.4 mm; 
a 0.5-arcmin source will be underestimated by 10% and 
11%; a 1-arcmin source by 31% and 36%. 



Given the ~ 1' beam of the SPT, we expect that 
few emissive sources will appear extended in our maps. 
With 1-arcmin resolution, a normal galaxy will appear 
point-like at redshifts z > 0.05 (distances greater than 
~ 200 Mpc), so only very nearby objects or objects with 
very extended structure (such as AGN with 100-kpc-scale 
jets) would appear extended in our maps. Furthermore, 
the matched filter applied to the maps is optimized for 
unresolved sources and will degrade the signal-to-noise 
on any extended source. 

We search for extended sources by fitting a cut-out of 
the (unfiltered) map around each detected source to a 
model of our measured beam convolved with a Gaussian 
of variable width. We then identify sources for which the 
best-fit FWHM is at least 0.25 arcmin and is inconsistent 
with zero at the 3a level. We also visually inspect the 
filtered map at each > 4.5ct source location for possible 
extended sources and any other anomalies. 

Of the 188 sources detected with S/N> 4.5 in cither 
band, 11 have a best-fit width of at least 0.25 arcmin and 
are inconsistent with zero width at > 3a. These sources 
are noted in Table 5 with an "a" next to the source 
name. Of these 11 sources, nine fall into our synchrotron- 
dominated class and have counterparts in the SUMSS or 
Parkes-MIT-NRAO (PMN, Wright et al. 1994) catalogs. 
Three of these nine are also listed in the SUMSS cata- 
log as having detectable extent beyond the ~ 30 arcscc 
SUMSS beam. The remaining two sources that we iden- 
tify as extended are nearby IR-luminous galaxies that 
are also detected with IRAS. Our visual inspection of all 
sources above 4.5cr in either band revealed the following 
cases of note (some of which are also identified by the 
quantitative test for extended structure): 

SPT-S J051614- 5429.6: This detection may be spu- 
rious, caused by sidelobes from the deep SZ decrement 
at 2.0 mm from the galaxy cluster SPT-CL J0516-5430 
(also RXCJ0516.6-5430 and Abell S0520). There is no 
counterpart at 1.4 mm or in external catalogs, the source 
is classified as extended by the method described above, 
and visual inspection shows it to have an irregular shape. 
The other bright source very near a galaxy cluster with a 
deep SZ decrement, SPT-S J050907-5339.2 (near SPT- 
CL J0509— 5342) is almost certainly not spurious, since 
it is detected more strongly at 1.4 mm (which is near 
the SZ null) than at 2.0 mm. Evidence of this source 
is also seen in data taken with the Atacama Cosmology 
Telescope (Hincks et al. 2009). 

SPT-S J051217-5724.0: This source is classified as 
extended, and visual inspection reveals a clear offset be- 
tween the 2.0 mm and 1.4 mm emission. The emis- 
sion in both bands is almost certainly associated with 
the low-redshift (z = 0.0047) galaxy NGC 1853; we ap- 
pear to be resolving different components of emission 
within the galaxy. We see a similar configuration in 
SPT-S J050656-5943.2, which is associated with the z = 
0.0041 galaxy NGC 1824, and SPT-S J055116-5334.4, 
which is associated with the z = 0.015 galaxy ESO 160- 
G 002. We also see offsets between 1.4 mm and 2.0 mm 
emission in SPT-S J051116— 5341.9 which has no obvious 
counterparts in existing catalogs. 

SPT-S j r 052850- 5300.3: We classify this source as 
dust-dominated, and it has no counterpart in the IRAS- 
FSC. There is a SUMSS source 45 arcsec away, and vi- 
sual inspection reveals a low-significance 2.0 mm counter- 
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Fig. 7. — Raw 1.4 mm flux versus raw 2.0 mm flux for sources detected above 4.5cr [dotted lines) in either band, long-dashed 
line: A spectral index a — 3 typical for sources dominated by dust emission, short-dashed line: A spectral index a = —1 typical 
for sources dominated by synchrotron emission. By finding associations within 1 arcmin for SUMSS (purple circles), AT20G 
(orange triangles), RASS (light blue squares), and IRAS (green diamonds), we see that most synchrotron-dominated sources are 
previously known. The bright dust-dominated population without counterparts in IRAS is discussed in Sec. 7.2. 
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part exactly coinciding with the SUMSS location. The 
1.4 mm emission, however, is clearly offset from both 
the 2.0 mm emission and the SUMSS location, indicat- 
ing that this may be a chance superposition of a known 
radio source and a previously unknown dust-dominated 
source. 

4. CORRECTING FOR FLUX BOOSTING AND 
ESTIMATING SPECTRAL BEHAVIOR 

The differential counts of mm-wave selected point 
sources as a function of source flux are expected to be 
very steep, so the measured flux of a point source in the 
SPT survey will almost certainly suffer flux boosting. In 
this work, we define flux boosting as the increased prob- 
ability that a source we measure to have flux S is really a 
dimmer source plus a positive noise fluctuation relative to 
the probability that it is a brighter source plus a negative 
noise fluctuation. Because of this asymmetric probability 
distribution, raw measurements of source flux will be bi- 
ased high. 29 The standard method in the SMG literature 
for dealing with this problem (e.g., Coppin et al. 2005) is 
to construct a posterior probability distribution for the 
intrinsic flux of each detection. The situation with SPT 
data is more complicated for two reasons: 1) As discussed 
in Crawford et al. (2010), the current implementation of 
this method in the SMG literature is not appropriate for 
estimating properties of individual sources, which is a 
key aim of this work; 2) We have data in more than one 
observing band, and the prior information that is ap- 
plied to create the posterior flux likelihood will be highly 
correlated in the two bands. 

In Crawford et al. (2010), we develop a method of cor- 
recting for flux boosting (based on the Bayesian poste- 
rior method used in Coppin et al. (2005) and others) 
which preserves information on individual source proper- 
ties, and we extend that method to estimate the intrinsic 
multi-band flux of a source based on the measured flux in 
each band and the prior knowledge of the source popu- 
lations in the various bands. In the two-band SPT case, 
the final product for each source is a two-dimensional 
posterior likelihood, where the two variables are either 
the flux in each band or the flux in one band and the 
spectral index between bands. The two likelihood distri- 
butions are trivially related by: 

)= (8) 
dct 

-P("Sniax,lj Ot\Sp t m,\i »Sp,m,2) > 

where S Pjm ,i is the measured flux in a resolution ele- 
ment or pixel in band i, SWiax,; is the true flux of the 
brightest source in that resolution element and band, and 
da/dS maXj 2 is derived from Eqn. 5. If we cast our prior 
information on source behavior in terms of source counts 
in one band and spectral behavior between bands, and 
we make the assumption that spectral index does not 

29 This phenomenon is closely related to what is referred to in 
the literature as "Eddington bias" (e.g., Teerikorpi 2004); however, 
the consensus use of the term in the literature is to describe the 
bias introduced to estimation of source counts vs. brightness, not 
on the estimated brightness of individual sources. This usage is 
consistent with the original work of Eddington (1913). As such, 
we choose to use "flux boosting" for the effect on individual source 
flux estimation. 



depend on flux, then we can write: 

P(Smax,l, &\Sp,m,l, Sp,m,2) (9) 
P{Sp, m ,l, <Sp,rn,2|<Smax,l, Ct)P(S nlaXt i)P(a) 

(see Crawford et al. (2010) for details). Of course, we 
expect that the spectral index distribution of our sources 
will in fact depend on flux. But the prior that we choose 
to place on a (see Sec. 4.1) is broad enough that it easily 
encompasses the full expected spectral index distribution 
at all fluxes. 

The posterior probability distributions in Eqn. 8 and 
Eqn. 9 are used to calculate most of the quantities re- 
ported in subsequent sections, including the 16%, 50%, 
and 84% (68% probability enclosed about the median) 
percentiles for de-boosted flux listed in Table 5; the prob- 
ability distributions for 2.0 mm-1.4 mm spectral index, 
from which the 16%, 50%, and 84% percentiles for that 
quantity in Table 5 are derived; and the source counts 
shown in Fig. 9, Fig. 10 and Fig. 12. 

4.1. Choice of Priors 

To construct the prior P(S maXt i) in Eqn. 9, we need 
assumptions for the source counts as a function of flux 
(dN/dS) in each of our bands. We assume the counts 
in each band will be the sum of dust-dominated and 
synchrotron-dominated counts, and we use estimates of 
counts for each of these populations from the literature, 
extrapolated to our wavelength bands when necessary. 
For the dust-dominated population, we use the Negrello 
et al. (2007) model for counts at 850 /jm, extrapolated 
to our wavelengths using a spectral index of 3.0 for the 
SMGs and 2.0 for the IRAS- type galaxies (assuming zero 
scatter in the index in both cases). The choice of these 
spectral indices was taken from an Arp220 SED tem- 
plate and the outcome is not very sensitive to the in- 
put. 30 De Zotti ct al. (2005) make direct predictions for 
the synchrotron-dominated population counts at 2.0 mm, 
which we use without modification. We extrapolate these 
predictions to 1.4 mm using a Gaussian distribution of 
spectral indices, centered on —0.5 with RMS of 0.5. We 
have found that the choice of source-count prior makes 
only a small difference in the resulting posterior proba- 
bility distributions (in the S/N range presented in this 
paper), consistent with the result in Scott et al. (2008). 
Sec. 2.6 of the companion paper Crawford et al. (2010) 
describes the interplay of experimental and prior infor- 
mation. 

For the spectral index prior, we have chosen a fiat prior 
between a = — 3 and a = 5. Given what is known about 
the two populations that are expected to contribute to 
sources at our wavelengths (e.g., Knox et al. 2004; Ma- 
son et al. 2009), this estimate conservatively brackets 
the expected spectral behavior of SPT sources. (The dif- 
ference in derived counts by choosing a flat prior from 
a = —2 to a = 4 is negligible.) Sec. 3.1 of the compan- 
ion paper Crawford et al. (2010) describes a subtlety of 
choosing the counts prior in parallel with a prior on the 
spectral index, because the two arc interdependent. In 

30 While we compare to the integral counts predictions of La- 
gache et al. (2004) in Sec. 7.2, we found that a kink in those dif- 
ferential counts produced a bias toward drawing ~ 10 mjy sources 
from the posterior. Negrello et al. (2007) was then used by virtue 
of its smoothness. 
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this work, we use a 1.4 mm counts model prior for the 
derived 1.4 mm counts and a 2.0 mm counts model prior 
for the derived 2.0 mm counts. This means that in the 
derived 1.4 mm counts, the information that comes from 
the 2.0 mm band is translated to 1.4 mm flux using the 
flat prior from a = -3 to a = 5 (and similarly for the 
cross-band information for the 2.0 mm counts). 



35 




Spectral index 



Fig. 8. — The distribution of the posterior spectral indices 
measured between 1.4 mm and 2.0 mm for sources with signal- 
to- noise > 4.5 in both bands (thick black line). Because we 
take a flat prior on the spectral index between —3 and 5 (and 
zero outside), the distribution outside the plotted range here 
goes to zero. We sum the P(a) for each source (described 
in Sec. 4.2 and normalize it such that the integral over a is 
equal to 1 for each source) to give an effective dN/ dot for this 
selection. We then classify the source by the probability that 
its posterior spectral index distribution exceeds a classifica- 
tion cut, taken here to be a = 1.66. (This threshold is at the 
minimum of dN/da between the two populations.) Sources 
with > 50% probability of posterior a > 1.66 are classified 
as dust dominated and those with < 50% probability of pos- 
terior a > 1.66 are synchrotron-dominated. There are 11 
dust sources (light gray) and 41 synchrotron sources (dark 
gray) that contribute to this distribution. The population 
split shown here is robust to changes in the signal-to-noise 
cut. At lower signal-to-noise cuts, the population features 
broaden slightly and many sources have poorly-localized P(a) 
distributions which contribute a floor in dN/da. 

4.2. Spectral Index Estimation and Source Classification 

Method 

By marginalizing the two-dimensional posterior in 
Eqn. 9 over the flux in the detection band P(S max ^i), 
we obtain a posterior likelihood for the spectral index of 
each detected source. The 16%, 50%, and 84% values 
of a given in Table 5 are taken from the cumulative ver- 
sion of this likelihood distribution for each source. These 
individual distributions can be summed to produce the 
measured a distribution of all sources detected in our 
two bands (which will be the convolution of the intrin- 
sic distribution with a complicated function of the noise 
from instrumental, atmospheric, and source background 
contributions in both bands). Fig. 8 shows that the pos- 
terior spectral index distribution for sources at S/N> 4.5 
in both bands has a clear population split. We use this 



split to identify the sources as either dust or synchrotron- 
dominated through the posterior. 

In Table 5, a source with P(a > 1.66) > 0.5 (hav- 
ing > 50% of its posterior index distribution in excess 
of 1.66) is classified as dust-dominated and a source 
with P(a > 1.66) < 0.5 (having < 50% of its poste- 
rior index distribution in excess of 1.66) is classified as 
synchrotron-dominated. The source counts by popula- 
tion use a probabilistic method based on P(a) that is 
described in Sec. 5, but there, too, we take a = 1.66 
to be the threshold. This value of the threshold is the 
spectral index at the minimum of dN/da between the 
two populations (see Fig. 8). The source counts pre- 
sented here are insensitive (within their uncertainty) to 
this choice of threshold over the range a = (1.66 ± 0.5); 
this particular value is chosen for definiteness. 

5. SINGLE-BAND SOURCE COUNTS 

Using the flux measurements in each band, we estimate 
the probability distribution of intrinsic fluxes for each 
source by constructing the two-band posterior likelihood 
P(S maK ,i, S maXt 2\S p ,m t i, S p , mt 2) as described in Sec. 4. 
To derive the counts as a function of flux from these 
distributions of intrinsic fluxes, we apply a bootstrap 
method similar to the one described in Austermann et al. 
(2010). Here, for each source, we randomly draw 5 x 10 4 
intrinsic fluxes from the two-band posterior, forming ef- 
fectively 5 x 10 4 mock catalogs of intrinsic fluxes in both 
bands. For each catalog, we draw a subset of these 
sources with replacement, where the number of sources 
drawn is a Poisson deviate of the catalog size. This re- 
sampling accounts for sample variance but not cosmolog- 
ical variance (which would require an additional variance 
term to describe how counts are expected to vary from 
sky patch to sky patch because of large scale structure) . 

For each of these 5 x 10 4 rcsamplcd catalogs, we esti- 
mate dN/dS and N(> S). In each flux bin, we then find 
the 16%, 50%, and 84% percentile points (that is, 68% of 
the enclosed probability around the median) of the dis- 
tributions of dN/dS and N(> S) in that bin. This yields 
the equivalent of la normally-distributed errors in each 
flux bin. The N(> S) and dN/dS bootstraps account 
for sample variance and posterior distribution variance 
(which includes noise, calibration error, deboosting, and 
cross-band information). Because the posterior flux dis- 
tributions per source span several flux bins in the counts, 
even the errors on differential counts will be correlated. 
The counts are corrected for the completeness using the 
simulations described in Sec. 3.3.1. Fig. 9 shows differ- 
ential counts in both bands as well as the differential 
counts of populations that we identify as either dust or 
synchrotron-dominated. The counts arc given in numer- 
ical form in Tables 1 through 4. 

The discrimination into synchrotron and dust- 
dominated counts is included in the bootstrap method. 
Here, each catalog resampling will have drawn intrinsic 
fluxes at 1.4 mm and 2.0 mm, and so each source in that 
catalog will have a spectral index. The dust-dominated 
(synchrotron-dominated) counts are derived from those 
sources that have a > 1.66 (a < 1.66) in the resampling. 
A source with P(a > 1.66) > 0.3 will therefore fall into 
the dust counts in 30% of the resamplings and into the 
synchrotron counts 70% of the resamplings. 

We also estimate purity in each flux bin of the counts 
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Fig. 9. — Differential source counts by population for the 2.0 mm (upper plot) and 1.4 mm (lower plot) bands. Gray boxes and 
black crosshairs indicate the total counts in that band. Red crosshairs indicate the synchrotron-dominated population counts 
and blue crosshairs indicate the dust-dominated counts. Crosshairs with full error enclose 68% of the probability about the 
median and are estimated in the bootstrap over flux described in Sec. 5. Here we have offset the two populations slightly in 
flux so that they do not lie on top of one another and the total counts are at zero offset. A source is identified as synchrotron- 
dominated (dust-dominated) if a < 1.66 (a > 1.66) in the bootstrap resampling from the joint posterior flux distributions, see 
Sec. 5. This splits the populations so that their differential counts sum to the total counts. A correction for survey completeness 
from simulations described in Sec. 3.3 is also applied, and impacts primarily the 1.4 mm counts in the lowest two flux bins. 
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statistics using the resamplcd catalogs. The purity is not 
evaluated for the flux at the bin center (which represents 
an intrinsic flux), but is instead related to the signal-to- 
noise of the raw flux of the sources that contribute to 
that bin. In each resampling, if a source lies in a flux 
bin, we find the associated purity from its raw signal-to- 
noise. The purity in the bin is then the weighted average 
of the purity of each source detection that contributes to 
the bin. In the counts presented here, the purity in the 
lowest flux bin is 0.80 at 2.0 mm and 0.76 at 1.4 mm. 

The Bayesian method accounts for a particular sense 
of the purity which is slightly different than purity pre- 
sented in Sec. 3.3.2. In Sec. 3.3.2, we take purity to be the 
fraction of noise fluctuations that are counted as sources. 
In principle, one would then suppress this fraction of the 
counts because it represents spurious detection. The out- 
look of the Bayesian method is that each pixel always 
has some intrinsic source flux - it may just be a tiny 
flux on top of a large positive noise fluctuation, and so 
is not strictly a false detection. When we determine the 
counts, we draw such sub-threshold fluxes from the pos- 
terior flux distribution. The effect of this is that in some 
fraction of the bootstrap samples, a source will scatter 
downward in flux (representing a tiny intrinsic flux plus 
a large positive noise fluctuation), out of the flux range 
presented. No additional correction for the purity (in the 
sense of Sec. 3.3.2) needs to be applied. Completeness 
is not included in this framework and so is corrected for 
explicitly. 

6. ASSOCIATIONS WITH EXTERNAL CATALOGS AND 
FOLLOW-UP OBSERVATIONS WITH ATCA 

Where possible, we identify candidate counterparts to 
the SPT-detected sources in several external catalogs and 
databases. We have queried the NED 31 and SIMBAD 32 
databases for counterparts within 2.0 arcmin of all 3496 
of our > 3cr sources. We have also searched catalogs 
from five individual observatories for counterparts: 1) 
the SUMSS catalog; 2) the IRAS-FSC; 3) the RASS-BSC 
and RASS-FSC; 4) the PMN catalog; and 5) the AT20G 
catalog. We search these catalogs, in particular, because 
these observatories are especially relevant for extragalac- 
tic sources in the Southern Hemisphere. Additionally, as 
mentioned in Section 3.3.2, we have performed follow- 
up observations at 6 cm on many of our brightest sources 
with the Australia Telescope Compact Array (ATCA). 

ATCA observations were performed at 6 cm, 12 mm, 
and 7 mm during August 2008 under program C1563. At 
6 cm, 55 sources in this field were observed and 52 sources 
were detected at a typical RMS noise of 0.5 mJy. One of 
the sources not detected (SPT-S ,1053250-5047.1) we clas- 
sify as a dusty source (see Section 7.2). The other two 
undetected sources (SPT-S J053412-5924.3 and SPT-S 
J055232-5349.4) are extended in the SUMMS catalog 
(45" resolution) and thus heavily resolved at the typi- 
cal 9" x 3" 6 cm resolution. 

The majority of our synchrotron-dominated sources 
have a clear counterpart either in an external catalog, 
in our 6 cm ATCA observations, or both. Because these 
sources are expected to have fluxes that vary significantly 
over time (Kellermann & Pauliny-Toth 1968), and be- 

31 http:/ /nedwww. ipac.caltech.edu/ 

32 http:/ /simbad.u-strasbg.fr/simbad 



cause the radio catalogs in which we search for coun- 
terparts are not significantly deeper than our own mea- 
surements, one might expect to not find counterparts for 
sources that were caught at peak brightness in the SPT 
observations. 

Despite this caveat, of the 107 sources above 5cr that 
we classify as synchrotron-dominated, only three of them 
do not have SUMSS counterparts within 30 arcsec. The 
brightest of these three sources (SPT-S J053345-5818.1) 
is classified as extended using the method of Sec. 3.5, 
and we find counterparts for this source within 10 arcsec 
in both the PMN catalog and our ATCA follow-up ob- 
servations. The extended nature of this source and the 
offset between the SUMSS and SPT/PMN/ATCA posi- 
tions for this source are expected if this source is an AGN 
with extended jets. The synchrotron emission from the 
source is presumably dominated by the radio-lobe (jet) 
contribution in the 36 cm SUMSS observations but dom- 
inated by emission from the core at shorter wavelengths. 
This frequency-dependent core-to-lobe flux ratio is com- 
monly seen in radio- loud AGN (e.g., Kharb et al. 2008; 
De Zotti et al. 2010) and is predicted by certain unified 
AGN models (e.g., Jackson & Wall 1999). Indeed, visual 
inspection of this source reveals one SUMSS source on 
either side of the SPT location (each within 40 arcsec) 
and a RASS-BSC object (also presumably dominated by 
emission from the AGN core) directly on top of the SPT 
location. 

Of the remaining two > 5a synchrotron-dominated 
sources with no SUMSS counterpart within 30 arcsec, 
one (SPT-S J050334-5244.8) has a SUMSS counterpart 
and a counterpart in our ATCA follow-up, both within 
35 arcsec, and the other (SPT-S J053726-5434.4) has a 
counterpart 11 arcsec away in the Veron-Cetty & Vcron 
(2006) and Hamburg-ESO (Wisotzki et al. 1991) catalogs 
as well as a possibly associated SUMSS source 1.3 arcmin 
away. We thus believe that every one of the sources that 
we detect at > 5c and classify as synchrotron-dominated 
is a real source. This is consistent with our estimates of 
purity in Sec. 3.3.2, which predict a false detection rate 
of effectively zero above 5a. 

The situation with our dust-dominated sources is very 
different. Of the 23 (47) sources above 5<7 (4.5c) that 
we classify as dust-dominated, only 10 (12) have coun- 
terparts (in any catalog) within 30 arcsec, and only 12 
(15) have counterparts within 1 arcmin. Given the stud- 
ies summarized in Sec. 3.3.2 and the counterparts found 
for the synchrotron-dominated sources, there is almost 
no chance that all (or even a majority of) these detec- 
tions without counterparts are spurious. Of the dust- 
dominated sources above 5c that do have counterparts 
within 30 arcsec, all but three are nearby galaxies de- 
tected with IRAS. Two of the remaining three are as- 
sociated with SUMSS sources, while the other is SPT- 
S J054716— 5104.1 which is associated with the debris 
disk around the star (3 Pictoris. In this field there are 
no other SPT sources within 30 arcsec of a SIMBAD 
database star. 

7. INDIVIDUAL-POPULATION SOURCE COUNTS AND 
IMPLICATIONS 

To briefly summarize the results of the last several sec- 
tions, three broad classes of point sources are detected 
with high significance in the SPT 1.4 mm and 2.0 mm 
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TABLE 1 
The 2.0 mm differential counts. 
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TABLE 2 
The 1.4 mm differential counts. 
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maps: 

1. Sources with 1.4 mm-to-2.0 mm flux ratios con- 
sistent with synchrotron emission, the vast ma- 
jority of which appear in radio catalogs and/or 
in our centimeter- wave follow-up observations with 
ATCA and which we generically refer to as AGN; 

2. Sources with 1.4 mm-to-2.0 mm flux ratios con- 
sistent with dust emission which have low-redshift 
(z -c 1) counterparts in the IRAS-FSC, and which 
we will generically call IRAS sources; 



3. Previously undetected sources with 1.4 mm-to- 
2.0 mm flux ratios consistent with dust emission. 

Here we present the measured source counts as a func- 
tion of flux for each of these populations and discuss the 
implications of our measurements. 

7.1. Synchrotron- Dominated Source Counts 

From associations with radio catalogs and from our 
6 cm followup observations with ATCA, we conclude 
that SPT sources with 1.4 mm-to-2.0 mm spectral in- 
dices less than 1.66 are consistent with being members of 
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TABLE 3 

The 2.0 mm cumulative counts. The reported purity is the weighted average of the 
purity of source detections that contribute to a bin, see sec. 5. 
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TABLE 4 
The 1.4 mm cumulative counts. 
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Fig. 10. — Differential counts for the population of sources identified as synchrotron-dominated compared to the De Zotti 
et al. (2005) model. Here the counts are scaled by S 2 ' 5 relative to by-band differential counts shown in Fig. 9 to match the 
(geometrical) convention in AGN literature. The error regions enclose 68% of the probability centered about the median counts, 
and are calculated using the bootstrap over the two-band posterior intrinsic flux (at 2.0 mm) that is described in Sec. 5. In 
these counts, a source is identified as synchrotron-dominated if a < 1.66 in the resampling. 
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the classical radio-source population (see Dc Zotti ct al. 
(2010) for a recent review). Although contributions to 
this population can come from synchrotron and free-free 
emission in normal and starburst galaxies (Condon 1992; 
De Zotti et al. 2010), the population is dominated at 
short radio wavelengths (30 cm and below) and moder- 
ate to high fluxes (10 mJy to 1 Jy) by synchrotron emis- 
sion from AGN (De Zotti et al. 2010). At even shorter 
wavelengths (1 cm and below), the moderate-to- high-flux 
counts arc expected to be dominated by the sub-class of 
AGN known as flat-spectrum radio quasars (FSRQs) (De 
Zotti et al. 2005). 

The behavior of this source population at mm wave- 
lengths is interesting for several reasons. Astrophysically, 
mm measurements of AGN have the potential to inform 
models of AGN emission mechanisms and evolution, par- 
ticularly whether FSRQs undergo spectral steepening at 
shorter wavelengths. This short-wavelength behavior of 
FSRQs is also of interest to the CMB and SZ commu- 
nities, as the predictions of contamination of mm CMB 
power spectrum measurements and SZ galaxy cluster sur- 
veys by AGN emission depend heavily on extrapolations 
of long- wavelength source properties to mm wavelengths 
(Lin & Mohr 2007; Rcichardt et al. 2009). Finally, the 
compact angular size of FSRQs (along with their short- 
wavelength brightness) make them attractive candidates 
for phase calibration sources for the Atacama Large Mil- 
limeter Array (ALMA). 

Two obvious questions that can be addressed with 
multi-band SPT observations of classical radio sources 
are: 1) Are the sources selected at 2.0 mm consistent with 
a population of FSRQs; and 2) Is the spectral behavior 
of these sources between 2.0 mm and 1.4 mm consistent 
with their longer- wavelength (flat-spectrum) behavior, or 
is there evidence of a spectral break or turnover? 

The first question is made somewhat ambiguous by the 
lack of a clear definition of the FSRQ population in terms 
of synchrotron spectral index. De Zotti et al. (2005) use 
a single spectral index of a = —0.1 for their model of the 
FSRQ population, while surveys targeting "flat spectrum 
radio sources" have generally chosen a > —0.4 or —0.5 as 
the defining threshold (Jackson et al. 2002; Healey et al. 
2007). Nevertheless, we can ask whether the behavior of 
these sources between 2.0 mm and longer wavelengths is 
more consistent with a mean spectral index such as De 
Zotti ct al. (2005) use or with a value more typical of 
steep-spectrum sources such as a ~ —0.8 (De Zotti et al. 
2010). 

We investigate this question using two lines of evi- 
dence. First, we have the comparison of the 2.0 mm 
SPT flux with fluxes measured at longer wavelengths. 
A serious caveat to any comparison of non-simultaneous 
observations of radio sources is that the known variabil- 
ity of such sources will add scatter — and, potentially, 
bias — to estimates of spectral behavior. Our ATCA 
6 cm follow-up observations took place approximately 
two months after the SPT observations, which will limit 
the effect of variability on our 6 cm-to-2.0 mm compar- 
isons to timescales shorter than this separation. Accord- 
ing to Murphy et al. (2010), the AT20G observations of 
sources between —60° and —50° declination took place 
in September and October, 2005, so 1.5 cm-to-2.0 mm 
comparisons should be viewed with greater caution. 

With this caveat in mind, we find that our highest- 



significance 2.0 mm-sclccted, synchrotron-dominated 
sources are completely consistent with flat spectral be- 
havior (or a ~ —0.1) between 6 cm and 2.0 mm. If we 
take the 57 sources for which we have robust (> 5a) de- 
tections at both 6 cm and 2.0 mm, and we calculate a 
single spectral index for each source using the best-fit 
2.0 mm flux from Table 5 and the raw 6 cm flux, we find 
these sources have a distribution of spectral indices char- 
acterized by a = — 0.13±0.21. Because this sample of 57 
sources is not purely 2.0 mm-selected (we are omitting 
sources without a robust 6 cm detection) and because the 
6 cm data have not been corrected for flux boosting, one 
might worry that this result is biased toward high 6 cm 
fluxes and, hence, steeper spectral indices. If we repeat 
the calculation using only the top 10 brightest 2.0 mm 
sources — which are all detected above 50tr at 6 cm — 
we find a spectral index distribution of a = —0.08 ±0.22. 

The second line of evidence that supports the hypoth- 
esis that our 2.0 mm-selected, synchrotron-dominated 
sources are FSRQs — or, more specifically, have a mean 
spectral index near a = —0.1 — is the agreement be- 
tween the 2.0 mm synchrotron-dominated counts and the 
predictions of De Zotti et al. (2005). The De Zotti et al. 
(2005) model includes contributions from many popula- 
tions of radio sources, including normal and star-forming 
galaxies and many types of AGN, but at 1 cm and below, 
the > 10 mJy model counts are dominated by FSRQs. 
Fig. 10 compares our synchrotron-dominated 2.0 mm 
counts to the De Zotti et al. (2005) 2 mm model. The 
model is roughly consistent with our measured counts, in- 
dicating that the 2.0 mm synchrotron-dominated counts 
are consistent with a power- law extrapolation of the long- 
wavelength FSRQ counts, and that the spectral behavior 
used to extrapolate the long- wavelength FSRQ counts in 
the model — in this case, simply assuming a = —0.1 for 
all FSRQs — is reasonably accurate down to 2.0 mm. 

The answer to the second question — whether the ra- 
dio sources selected at 2.0 mm show flat spectral be- 
havior all the way to 1.4 mm — is addressed by our 
simultaneous 1.4 mm and 2.0 mm observations of these 
sources. These observations indicate that a mean spec- 
tral index of a = —0.1 is not an accurate description of 
the spectral behavior of these sources between 1.4 mm 
and 2.0 mm. The distribution of spectral indices in 
Fig. 8 shows that our synchrotron-dominated sources 
have 1.4 mm-to-2.0 mm spectral indices peaked around 
a = —0.5, and a similar result is found for the sum of the 
posterior spectral index PDFs for the 57 sources detected 
above 5a at both 6 cm and 2.0 mm. If we take the ten 
brightest 2.0 mm-selected sources, we find that the sum 
of the posterior spectral index PDFs of these sources is 
inconsistent with a = —0.1 at 97%. 

Interestingly, if we perform the same calculation us- 
ing the published 1.5 cm AT20G fluxes for our top ten 
2.0 mm-selected sources, the 2.0 mm-to-1.5 cm spectral 
index distribution is a = —0.31 ± 0.29 while the 1.5 cm- 
to-6 cm spectral index distribution is a = 0.25 ± 0.36 — 
i.e., there is some evidence that the flux of these sources 
peaks between 6 cm and 2.0 mm. Biases induced by 
source variability are unlikely to have a significant effect 
on this measurement for two reasons. First, these sources 
are so far above the detection threshold at 2.0 mm that 
they would have been significant detections at any epoch. 
Second, biases due to variability would tend to drive 
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estimates of this spectral index in the opposite direc- 
tion — because sources are likely to be selected at their 
peak brightness, measurements taken at any other epoch 
should produce systematically lower fluxes. This implies 
that if we had 1.5 cm measurements taken simultane- 
ously with the SPT measurements, we would tend to 
find even higher 1.5 cm fluxes than in the AT20G cata- 
log (assuming the variability is simultaneous at the two 
wavelengths). These arguments imply that the evidence 
for these sources' flux peaking between 6 cm and 2.0 mm 
is fairly robust, if not overly statistically significant. 

7.2. Dust- Dominated Source Counts 

7.2.1. The Dust-Dominated Assumption, Spectral Indices, 
and Self-Absorbed Synchrotron 

We have referred throughout this work to the popula- 
tion of sources with 1.4 mm-to-2.0 mm spectral indices 
greater than 1.66 as "dust-dominated." For the sources in 
this class that have IRAS-FSC counterparts, the assump- 
tion that thermal emission from dust is the dominant 
mechanism at mm wavelengths is reasonable given their 
large IR fluxes. The 1.4 mm-to-2.0 mm spectral indices 
for these sources are consistent with dust emission — the 
sum of the posterior spectral index PDFs for the 5 bright- 
est sources with IRAS counterparts peaks at a — 3.2 — 
and comparisons of 1.4 mm SPT fluxes and 100 IRAS 
fluxes of these sources show that their emission is con- 
sistent with thermal dust at moderate temperatures (20 
to 40 K) from the mm through the far-IR (see Sec. 7.2.2 
and Fig. 11 for details). 

While thermal dust emission is also the most natu- 
ral candidate for explaining the sources with 1.4 mm- 
to-2.0 mm spectral indices greater than 1.66 that do 
not have IRAS-FSC counterparts, we cannot a priori 
rule out self-absorbed synchrotron emission from AGN 
as the dominant emission mechanism for these sources. 
Self-absorbed synchrotron is the leading emission model 
for the population of gigahertz-pcaked-spectrum (GPS) 
radio sources (see O'Dea 1998 for a review), and it is 
the only emission mechanism other than thermal emis- 
sion from dust that could plausibly produce 1.4 mm-to- 
2.0 mm spectral indices well above 0. As the name sug- 
gests, GPS sources typically peak around wavelengths of 
30 cm (e.g., Stanghellini et al. 1998), but GPS sources 
have been observed with peaks at wavelengths as short 
as 1 cm (Edge et al. 1998), and there is no fundamental 
physics that rules out self-absorbed synchrotron emission 
peaking at much shorter wavelengths. However, several 
lines of reasoning suggest that the a > 1.66 sources with- 
out IRAS-FSC counterparts are not dominated by self- 
absorbed synchrotron. 

The first argument against GPS radio sources as the 
explanation for our "dust-dominated" counts without 
IRAS-FSC counterparts is that their spectral behavior 
is too steep even for self-absorbed synchrotron, while 
it is perfectly consistent with thermal dust emission. 
Stanghellini et al. (1998) show that, even well longward 
of the peak wavelength, the mean spectral index of GPS 
sources is a ~ 0.8, and rarely do they find spectral in- 
dices as high as 2.0. In contrast, the sum of the posterior 
spectral index PDFs for the brightest 5 dust-dominated 
sources without IRAS counterparts peaks at a = 3.3 and 
is inconsistent with a = 2.0 at 97% confidence. This peak 
value of a = 3.3 is consistent (within the width of the 



two distributions) with the peak of a = 3.2 for the top 
five sources with IRAS-FSC counterparts, sources which 
we are confident are truly dust-dominated. The peak 
value of these a distributions are also consistent with 
the best-fit mean spectral index of a = 3.7 ± 0.2 found 
by Hall et al. (2010) for Poisson-distributed sources be- 
low the SPT detection threshold — the vast majority 
of which are expected to be dust-dominated. Interest- 
ingly, all of these peak a values are actually steeper than 
typical values predicted by models for spectral indices 
of dust-dominated sources (e.g., Lagache et al. 2004), a 
point discussed in detail in Hall et al. (2010). 

Another argument against the GPS explanation for 
these sources is the lack of radio and x-ray counterparts. 
Sicmiginowska et al. (2008) found that GPS sources have 
2-10 keV fluxes of up to 10 46 erg s _1 , easily detectable 
in the ROSAT All-Sky Survey, and our brightest "dust- 
dominated" sources would have to be almost an order of 
magnitude dimmer at 36 cm than at 1.4 mm to evade 
detection in SUMSS. Finally, Kellermann & Pauliny- 
Toth (1981) argue that the peak wavelength of a GPS 
source should be proportional to flux density to the —0.4 
power, meaning that sources that peak at mm wave- 
lengths should be 2.5 orders of magnitude brighter than 
sources that peak at cm wavelengths, so they should be 
much rarer as well. Based on this set of arguments, we 
conclude that our sources with 1.4 mm-to-2.0 mm spec- 
tral indices greater than 1.66 are indeed dominated by 
thermal dust emission. 

7.2.2. SPT Dust-Dominated Source Counts and Arguments 
for a New Population of mm Sources 

One might presume that the dust-dominated sources 
presented here are exact analogues to the submm- 
selected SMG population measured by SCUBA and other 
instruments. However, SMG source counts at 850 /jm are 
found to drop precipitously at flux levels above 5-10 mJy 
(e.g., Coppin et al. 2006). Assuming an average spectral 
index of a = 2.8 (the estimate for mean SMG spectral 
index from Knox et al. (2004)), this drop in the counts 
would occur at ~ 2 mJy at 1.4 mm, leading one to ex- 
pect a 1.4 mm survey to see very few SMGs above the 
~ 10 mJy level needed for a robust detection in the SPT 
1.4 mm band. Indeed, extrapolating two of the mod- 
els used to fit SCUBA counts in Coppin et al. (2006) to 
1.4 mm and 10 mJy indicates that there should be <C 1 
SMG above 10 mJy at 1.4 mm in the SPT 87 deg 2 field; 
a similar extrapolation of the Schechter function fit to 
the AzTEC 1.1 mm counts in Austermann et al. (2010) 
yields a similar prediction. Using a spectral index closer 
to the mean of the dust-dominated part of the distribu- 
tion shown in Fig. 8 would drive these predictions even 
lower. Quite contrary to these predictions, SPT detects 
20 dust-dominated sources above 10 mJy at 1.4 mm, in- 
cluding 9 above 20 mJy. What are these sources if not 
SMGs? 

One possibility is that they are nearby ULIRGs, the 
low-redshift analogues of SMGs. These sources are rare 
enough to not contribute significantly to the submm 
counts estimated from very small patches of sky. As 
noted in Sec. 6, a fraction of the bright dust-dominated 
sources detected by the SPT have counterparts in the 
IRAS-FSC, and most of these sources are indeed low- 
redshift ULIRGs. The majority of SPT dust-dominated 
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Fig. 11.— IRAS 100 /an flux vs. SPT 1.4 mm flux for all SPT dust-dominated sources at S/N > 5. IRAS flux is taken 
from a version of the ISSA (Wheelock et al. 1994) 100 /im map which has been filtered to enhance point-source signal-to-noise. 
Horizontal error bars are the 68% enclosed interval in the posterior 1.4 mm flux distribution (as in Table 5). Vertical error bars 
are the width of the noise distribution in the filtered IRAS map. SPT sources with counterparts within 1 arcmin in the IRAS 
FSC are shown with diamond symbols. Lines of constant 100 [im-lA mm flux ratio are shown for five emission models, all 
modified blackbody laws with a dust emissivity index of /3 = 1.5 (consistent with the value of /3 used in Dunne & Eales (2001), 
Chapman et al. (2005), and Kovacs et al. (2006)) and with dust temperatures of 15, 20, 25, 30, and 35 K (if the emitter is nearby) 
or those temperatures times 1 + z (if the emitter lies at redshift z). There is a clear distinction between the locus of sources with 
IRAS-FSC counterparts — which have flux ratios consistent with warm, low-redshift dust emission — and the points which lie 
along the x axis and have no counterparts, which must be either at moderate to high redshift or have anomalously cold dust. 
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sources, however, do not have IRAS-FSC counterparts, 
or counterparts in any existing catalog. Likewise, in deep 
(~ 24.4 AB mag) griz optical data taken by the Blanco 
Cosmology Survey (BCS) 33 (which encompass roughly 
half of the 87 deg 2 described here) there are no obvious 
counterparts. 

It is unlikely that these sources are just below the 
threshold for inclusion in IRAS-FSC, given that the 
brightest three SPT dust-dominated sources do not have 
IRAS-FSC counterparts. Fig. 11 makes an even stronger 
argument that these sources are a very different popula- 
tion than the SPT sources with IRAS-FSC counterparts. 
This figure shows a scatter plot of IRAS 100 /xm flux — 
estimated from an IRAS Sky Survey Atlas (ISSA, Whee- 
lock et al. 1994) 100 fim map filtered to enhance point- 
source S/N — vs. SPT 1.4 mm flux. 34 The sources with 
and without IRAS-FSC counterparts occupy clearly dis- 
tinct loci of points. The sources with IRAS-FSC counter- 
parts are consistent with nearby sources that have typical 
dust temperatures of 25 — 35 K (Dunne & Eales 2001), 
while the sources without IRAS-FSC counterparts have 
ISSA 100 fim flux consistent with zero and thus must ei- 
ther be at moderate to high rcdshift or have anomalously 
cold dust (< 15 K). A stacking analysis of the 13 objects 
classified as dust-dominated with 1.4 mm S/N greater 
than 4.5 and no IRAS-FSC counterpart shows that these 
sources have mean IRAS 100 fim flux of 15.9 mJy, with 
a 2(7 upper limit for that mean of 106 mJy. These same 
sources have mean 1.4 mm flux of 18 mJy. If we assume 
an Arp220 SED for these sources, this stacking analysis 
implies these sources lie at a mean redshift of (z) =3.1, 
with a 2a lower limit of (z) = 2.0. Assuming a modi- 
fied blackbody spectrum with /? = 1.5, we obtain a 2a 
upper limit of the combination of dust temperature and 
redshift of Td/ (1 + z) = 14. 3K, implying a 2a lower limit 
on mean redshift of (z) = 1.1 if Td = 30K. 

The mean dust temperature of high redshift SMGs 
measured from previous work is ~ 35 K (Chapman et al. 
2005; Kovacs et al. 2006), but there have been detec- 
tions of a few galaxies with dust temperatures ~ 20K 
(e.g. Kovacs et al. (2006)). There is expected to be 
both a hot (~ 40K) dust component from the ISM sur- 
rounding actively star forming regions with many young 
stars, as well as a cold (~ 20K) dust component sur- 
rounding the diffuse quiescently evolving population of 
old red stars (Dunne & Eales 2001; Vlahakis et al. 2005; 
Coppin et al. 2008). IRAS would be largely insensitive 
to the cold dust component, which might explain why 
these SPT sources escaped previous detection. However, 
given the brightness of these objects, if these sources were 
galaxies at z < 1 then we would expect to see them in 
BCS, Digitized Sky Survey, 35 or Two Micron All Sky 
Survey images, 36 and we do not. This test depends on 
the assumption of a particular SED for the galaxies, but 
the most pessimistic galaxy SED assumption — namely 
an Arp220-likc SED with a high mm-to-optical bright- 
ness ratio — still predicts that such objects would be 

33 http:/ /cosmology.uiuc.cdu/BCS/ 

34 We chose to plot the IRAS 100 fim channel (as opposed to 
60 fim) because it is the closest band to the SED peak and to the 
SPT 1.4 mm band. 

35 http:/ /cadcwww.dao.nrc.ca/cadcbin/getdss 

36 http:/ /irsa. ipac.caltech.edu/applications/2MASS/ 



detected in the optical and near-IR (NIR). 

Throughout this work we have assumed that the ma- 
jority of sources detected in the mm waveband by the 
SPT have been of extragalactic origin, as this field is at 
high galactic latitude and has explicitly been chosen for 
its low contamination by galactic cirrus. How robust of 
a statement can be made as to the extragalactic nature 
on the dusty sources not detected in any external cata- 
log? A similar question was posed in Lawrence (2001) 
and it was demonstrated that the majority of SCUBA 
sources could not be of galactic origin. The SPT sources 
cannot be galactic HII regions as the accompanying free- 
free emission would have been detected by low frequency 
radio catalogs such as PMN and SUMSS. They could 
be nearby, cold (Td < 10 A), pre-star-forming clouds in 
the galaxy similar to sources described in Desert et al. 
(2008). However, as this field is at high galactic lati- 
tude and the clouds would presumably be located in the 
galactic disk and not the halo, the sources would there- 
fore be nearby (within ~ 1 kpc). Since these clouds are 
expected to be parsec-scale objects (Egan et al. 1998), 
they would have an angular extent (~ 5') much larger 
than the SPT beam and would be obviously detected 
by eye in the unfiltered map. We have also tested the 
dark cloud hypothesis by measuring optical extinction 
along the lines of sight towards the SPT dusty sources. 
A comparison of the optical color distribution of BCS 
objects toward dust-dominated sources to the colors to- 
ward random directions in the field shows no excess ex- 
tinction at the positions of these sources. Although we 
cannot definitively confirm that none of these sources are 
of galactic origin, these tests demonstrate that they are 
predominantly extragalactic. 

Comparisons of SPT dust-dominated source counts 
with model predictions strengthen the conclusion that 
our bright, dust-dominated sources without IRAS coun- 
terparts are a new, possibly high-rcdshift and lensed, 
population. Fig. 12 shows the cumulative SPT dust- 
dominated counts vs. flux and predictions from three 
models. Lagache et al. (2004), Negrello et al. (2007), 
and Pearson & Khan (2009) make predictions for counts 
at or very near our 1.4 mm band. All three models 
have two basic components: modcratc-to-high-rcdshift 
starburst galaxies (which account for basically all the 
counts seen by SCUBA at 850 fj,m) and nearby galax- 
ies (including the LIRGs and ULIRGs seen in IRAS). 
The Negrello et al. (2007) model also includes a com- 
ponent of the high-rcdshift starburst population which 
has been strongly lensed by foreground galaxies. At 
first glance, all models agree fairly well with the SPT 
counts at 1.4 mm. However, the counts at fluxes above 
10-20 mJy arc dominated in all models by sources that 
should be detectable in the IRAS-FSC above the 60 fim 
flux cut of 200 mJy, while our measured counts are dom- 
inated by sources without IRAS-FSC counterparts. We 
have modified the publicly available Lagache et al. (2004) 
code to exclude such sources from their model, 3 and 
Pearson & Khan (2009) and Negrello et al. (2007) have 
supplied us with model counts excluding sources with 

37 Using the template SED models supplied by Lagache et al. 
(2004), we find that any source detected at > 10 mJy at 1.4 mm 
below redshift z = 0.2 should have 60 (im flux above 200 mJy, so 
our modification to their model is effectively just a rcdshift cut at 
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Fig. 12. — Cumulative SPT dust-dominated source counts vs. 1.4 mm flux, with models overplotted. Model curves are as 
follows: short dashed line: Lagache et al. (2004) 1.38 mm prediction; long dashed line: Negrello et al. (2007) 1.4 mm prediction; 
dot-dashed line: Pearson & Khan (2009) 1.38 mm prediction. The error regions enclose 68% of the probability centered about the 
median counts, and are calculated using the bootstrap over the two-band posterior intrinsic flux (at 1.4 mm) that is described in 
Sec. 5. In these counts, a source is identified as dust-dominated if a > 1.66 in the resampling. Top Panel: This plot shows counts 
and models with all dust-dominated sources included. Bottom Panel: This plot shows counts calculated excluding sources that 
have IRAS-FSC counterparts (within 1 arcmin) and models calculated excluding populations that should be detectable in the 
IRAS 60 (im band above the typical FSC limit of 0.2 Jy. The SPT detects sources in excess of what is predicted by the models 
shown once the IRAS sources have been removed. 
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60 /jm flux greater than 200 mJy. We then re-calculate 
our 1.4 mm dust-dominated counts excluding sources 
with IRAS-FSC counterparts and compare these mod- 
ified counts to the modified predictions in the bottom 
panel of Fig. 12. 

There are significant discrepancies between our mea- 
sured counts without IRAS counterparts and the La- 
gache et al. (2004) and Pearson & Khan (2009) model 
predictions in both the slope of the counts and the num- 
ber of sources above 20 mJy. The Negrello et al. (2007) 
model also under-predicts the SPT counts but the shape 
of the counts are in good agreement. A modest adjust- 
ment along either axis would bring the model into ex- 
cellent agreement. As stated above, the main difference 
between the Negrello et al. (2007) model and the Lagache 
et al. (2004) and Pearson & Khan (2009) models is the 
inclusion of a population of high redshift, strongly lensed 
dusty star forming galaxies. 

Among the possible explanations for this new popula- 
tion, the hypothesis that these sources are at high red- 
shift is particularly compelling. The longer-wavelength 
SPT observations will be sensitive to a higher-redshift 
population than the submm surveys (due to the stronger 
negative K-correction), and there is considerable evi- 
dence that the very brightest SMGs are at the highest 
redshifts (Ivison et al. 2002; Pope et al. 2005; Greve et al. 
2008). This empirically observed trend of SMG bright- 
ness with redshift is plausible both because more distant 
systems have a higher probability of being gravitationally 
lensed (Blain 1996; Blain et al. 1999) and because evolu- 
tion in star formation as a function of environment, called 
"cosmic downsizing" , is consistent with a higher star- 
formation rate in massive systems at high redshift (Cowic 
et al. 1996; Cole et al. 2000). Models exist in the liter- 
ature that predict the existence of a mm/submm-bright 
population similar to ours resulting from both lensing 
(Blain 1996; Perrotta et al. 2003; Negrello et al. 2007) and 
intrinsically luminous high-z sources (Devriendt et al. 
2010). The reason that such objects would have been 
missed by previous mm/submm instruments is simply 
that SPT surveys so much more area (87 square degrees 
for this small subset of the SPT survey to ~ 1 square 
degree for the total area surveyed by SCUBA) and is 
hence much more likely to find rare, bright systems (due 
to strong lensing or intrinsic luminosity) that a smaller 
survey might miss. 

If this subset of the SPT-identificd dusty sources are 
indeed at high redshift, they represent an intriguing new 
class of mm sources, whether they are strongly lensed 
or intrinsically ultra-luminous. Strongly lensed systems 
allow observers to detect fainter background sources at 
higher redshift than would otherwise be obtainable. Be- 
cause lensing is achromatic, these sources will be brighter 
at all wavelengths, facilitating detailed studies which 
have otherwise been difficult to achieve. Equally excit- 
ing would be the identification of high-redshift galaxies 
which arc more massive and forming stars more prodi- 
giously than any systems yet identified. The identifica- 
tion of such galaxies would be a strong test of models of 
galaxy formation and evolution. Regardless of whether 
the sources are lensed or intrinsically luminous, a sample 
of high- redshift, dust-enshrouded star- forming galaxies 
has the potential to be a useful tool for the study of very 
early epochs of star and galaxy formation. 



8. CONCLUSIONS 

The SPT has detected 188 sources above 4.5c (over 
3000 above 3cr) in two-band data over a small (87 deg 2 ) 
subset of the full survey region. This is the first survey 
of its kind at mm wavelengths and mJy flux levels, both 
in survey area — the 87 deg 2 presented here is over an 
order of magnitude more area than previous mm surveys 
at these flux levels — and in the ability to spectrally dis- 
criminate between different source populations. Some of 
the sources detected in this work appear to be members 
of a new and intriguing population of dust-enshrouded, 
star-forming galaxies. 

We use the ratio of flux in the 1.4 mm and 2.0 mm 
bands to estimate the spectral index for every detected 
source. Using the posterior PDF for this spectral 
index, we classify each source as either synchrotron- 
dominated or dust-dominated. At high flux levels (above 
~ 15 mJy) in both bands, the majority of sources 
we detect are synchrotron-dominated; the synchrotron- 
dominated population continues to dominate the 2.0 mm 
counts down to the detection threshold, but the dust- 
dominated sources begin to take over the 1.4 mm counts 
near the 15 mJy level. 

The synchrotron-dominated sources we detect are con- 
sistent with the population of radio sources that is well- 
established from radio and cm surveys. All of our > 5er 
sources in this class have a clear counterpart in exist- 
ing radio catalogs. The number counts as a function of 
2.0 mm flux that we derive for this population are con- 
sistent with the predictions of De Zotti et al. (2005). 
Spectral comparisons within the SPT data and with 
ATCA data at 6 cm (taken in follow-up observations) and 
1.5 cm (Murphy et al. 2010) show that the synchrotron- 
dominated sources we detect are consistent with FSRQs 
(as predicted), but that the flat-spectrum behavior has 
turned over by mm wavelengths. This conclusion is im- 
portant for AGN models and for predictions of radio 
source contamination in CMB and SZ measurements. 

A fraction of our dust-dominated sources have coun- 
terparts in the IRAS-FSC and are typically associated 
with low-redshift (z -C 1) ULIRGs. The majority of our 
dust-dominated sources, however, have no counterpart in 
existing catalogs, and we argue that they represent a new 
and exciting population of mm sources. Comparisons of 
source counts for this population with model predictions 
and IRAS limits on the 100 /im flux of these sources 
demonstrate that these sources are inconsistent with ei- 
ther simple extrapolations of the submm-selected popu- 
lation at lower fluxes or with low-redshift galaxies with 
normal dust properties. Further, the mm-wavelength se- 
lection of such high-flux sources over a large survey area 
suggests a high-redshift population. Possibilities for ex- 
plaining this new family of sources include strong lensing 
of dimmer background sources, ultra-luminous starburst 
galaxies at moderate-to-high redshift, and (z ~ 1) galax- 
ies with extremely cold dust. Comparisons to models 
favors the hypothesis for these sources being strongly 
lensed. Any one of these explanations would have in- 
teresting implications for models of star and galaxy for- 
mation and be of potential cosmological interest. Prop- 
erly locating these sources in the broader context of the 
IR/submm/mm galaxy population is a major challenge 
for future work. An extensive program to determine the 



2G 



Vieira, et al. 



SEDs, morphologies, and rcdshift distribution of these 
objects is underway. Multi-wavelength follow-up imag- 
ing facilities such as the Hubble Space Telescope, the 
Spitzer Space Telescope, the Hcrschcl Space Observatory, 
and the Atacama Large Millimeter Array, will be crucial 
to disentangling the various possibilities and uncovering 
the nature of these objects. 

The point source results presented here use only a small 
fraction of the complete SPT data set. SPT is continuing 
to take data and has already observed 800 deg 2 to similar 
depths as this work at 1.4 and 2.0 mm. 600 deg 2 of 
this area also has 3.2 mm coverage. The complete SPT 
survey is expected to cover over 2000 deg 2 with these 
three wavelengths and will produce a catalog containing 
thousands of additional sources. 
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TABLE 5 

Point sources above 4.5ct at 1.4 mm or 2.0 mm in an 87 square degree SPT field 

CENTERED AT R.A. 5 h 30"\ DECL. -55° (J2000) 



ID &: coordinates: 2.0 mm data: 1.4 mm data: Spectral index & type: Nearest source in: 

SPT ID RA DEC S/N S raw S dist S/N S raw S dist cv raw ct dist P(q > 1.66) Type SUMSS RASS IRAS 

[deg] [deg] [ m Jy] [ m Jy] [ m Jy] [ m Jy] [arcscc] [arcscc] [arcscc] 



SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 



J045913- 
J050000- 
J050003- 
J050019- 
J050211- 
J050321- 
J050329- 
J050334- 
J050401- 
J050424- 
J050437- 
J050508- 
J050511- 
J050523- 
J050526- 
J050528- 
J050607- 
J050617- 
J050620- 
J050624- 
J050656- 
J050658- 
J050716- 
J050725- 
J050732- 
J050747- 
J050758- 
J050831- 
J050847- 
J050907- 
J051003- 
J051018- 
J051029- 



5942.4 

5752.5 

5229.0 

5321.3 

5040.7 

5328.3 

5735. 6 a 

5244. 8 b 

5023.2 

5711.9 

5818.1 

5346.4 

5346.0 

5808.5 

5044.8 

5056.6 

5844.8 

5748.7 

5741.0 

5024.8 

5943.2° 

5435.0 

5954.8 

5013.3 

5104.2 

5156.1 

5850. 5 a 

5449.4 

5754.3 

5339. 2 d 

5651.9 

5719.7 

5624.8 



74.804 
75.000 
75.015 
75.081 
75.546 
75.839 
75.874 
75.892 
76.006 
76.100 
76.156 
76.284 
76.297 
76.346 
76.358 
76.369 
76.533 
76.573 
76.584 
76.600 
76.735 
76.742 
76.820 
76.855 
76.884 
76.948 
76.995 
77.132 
77.199 
77.283 
77.516 
77.578 
77.623 



-59.708 
-57.875 
-52.485 
-53.356 
-50.679 
-53.472 
-57.595 
-52.747 
-50.387 
-57.199 
-58.302 
-53.775 
-53.767 
-58.143 
-50.747 
-50.944 
-58.748 
-57.813 
-57.685 
-50.415 
-59.721 
-54.585 
-59.914 
-50.223 
-51.071 
-51.936 
-58.842 
-54.823 
-57.907 
-53.653 
-56.865 
-57.329 
-56.414 



6.20 
44.01 
5.33 
39.42 
6.76 
5.56 
45.35 
5.96 
40.49 
16.05 
4.83 
4.76 
2.06 
4.56 
12.82 

I. 71 
7.25 

21.73 
5.61 
5.37 
5.18 
5.15 
4.91 
7.16 

19.68 

II. 75 
47.34 
0.84 
0.71 
3.27 
10.93 
23.24 
6.37 



7.51 
53.28 
6.70 
49.52 
8.86 
6.98 
54.90 
7.49 
53.07 
20.17 
5.85 
5.98 
2.59 
5.52 
16.81 
2.24 
8.77 
26.30 
6.79 
7.04 
6.27 
6.46 
5.94 
9.38 
25.79 
15.40 
57.31 
1.06 
0.86 
4.11 
13.74 
29.19 
8.00 



7 26+ 1 ' 27 
52.99±=;JS 

5.9ill;l? 

49.26±=;g 
8.16±i;£ 
6-1611:11 

54.60j^'j£ 
6.80±lj| 

52.79± 2 /.f 9 

19.94+i;« 

4.66+J;£ 

9 79 + 1.08 
£ - '^-1.42 

1 44+ 3 ' 44 

1.19 
16.507.! .48 

i.63ll;tf 

• 28 -1.30 
26.06±i;|l 

R+1.31 

-1.35 
c + 1.44 
-1.56 
A + 1.25 
-1.29 
0+1.37 
-1.42 

rj+1.53 

-3.02 

n+1.39 
-1.41 

25.56^! 69 

R +1.47 
-1.44 
! +2.58 
-2.58 

r,+ 1.08 

0.50 
+ 0.96 
0.44 
+ 1.30 
-1.33 

I3.4iti;!» 

+ 1.75 
1.70 

7.30t 1 1 ;f 6 



6.16j 
6.16j 
5.94j 
5.72j 
4.62 4 
8.70j 
25.56 
15.06 
57.01 
0.72 
0.64 
3. 



28.93± 



7.05 
15.15 
2.15 
10.82 
1.15 
1.43 
16.43 
1.95 
12.50 
7.56 
0.81 
1.49 
5.39 
0.22 
4.81 
4.87 
2.54 
6.01 
3.25 
1.64 
5.22 
3.64 
1.10 
1.09 
7.39 
3.46 
14.92 
4.66 
4.56 
4.61 
3.97 
5.68 
1.28 



22.58 
48.52 

7.20 
36.21 

3.98 

4.77 
52.62 

6.51 
43.27 
25.28 

2.59 

4.98 
18.05 

0.72 
16.65 
16.84 

8.12 
19.26 
10.42 

5.68 
16.73 
12.18 

3.53 

3.76 
25.59 
11.97 
47.79 
15.58 
14.60 
15.42 
13.30 
19.01 

4.27 



20.85lJ;»S 
47.80tlf a 

+ 2.88 
2.47 
+4.64 
4.64 
+ 2.51 
1.48 
+ 2.64 
1.70 
+ 5.76 
5.66 
+ 2.86 
2.19 
+ 5.17 
5.27 



5.57 
35.52 

4.71 

4.35 
51.83 

5.43 
42.54 



24 r;o+4.08 



2.69 
3.54 



+ 2.20 
1.66 
+ 2.70 
-1.95 



10 64 + 575 



1.57 
15.65 
1.94 
6.93 
18.47 
9.00 
4.73 
13.63 
10.79 
3.14 
4.79 
24.79 
10.96 
47.01 
1.10 
1.03 
4.76 
12.19 
18.20 
4.52 



+ 1.86 
1.26 



+ 3.78 
3.84 
+ 9.28 
1.29 
+ 3.03 
2.65 
+ 3.70 
3.65 
+ 3.26 
3.42 
+ 2.81 
1.98 
+4.17 
6.43 
+ 3.50 
-3.70 
+ 2.40 
1.57 
+ 2.43 
1.41 
+4.21 
4.16 
+ 3.64 
3.32 
+ 5.37 
i.27 
+ 2.02 
0.67 
+ 1.75 
0.63 
+ 9.05 
3.24 
+ 3.55 
3.55 
+3.78 
3.72 
+ 2.55 
1.54 



3.0 


_ __|_n 7 
^ • y - - 7 


0.95 


dust 


43 


575 


815 


-0.3 


-0.3t°oi 


0.00 


sync 


3 


313 


486 


0.2 


^ ^ + 1 4 


0.08 


sync 


18 


923 


1979 


-0.9 


n r^A-Cl 3 

-0.9±oI 


0.00 


sync 


4 


627 


189 


-2.2 


, - + 1 2 

-l-5-i.o 


0.00 


sync 


7 


1184 


1459 


— 1.0 


-0.9±J;1 


0.03 


sync 


8 


694 


1312 


-0.1 


,-, i -1-0 3 

— n 1 ' 


0.00 


sync 


2 


1058 


583 


-0.4 


^ - + 1 3 

-0-5-1.1 


0.02 


sync 


32 


581 


1611 


-0.6 


-0-6±g;i 


0.00 


sync 


1 


998 


1301 


0.6 


U '°-0.5 


0.00 


sync 


11 


158 


493 


-2.2 


^ ,-,-1-2 1 


0.13 


sync 


36 


895 


211 


-0.5 


„ „-|-2 


0.15 


sync 


120 


744 


605 


5.3 


n ^+1 n 
3.6_„ 


0.74 


dust 


82 


766 


508 


-5.6 


„ ^ + 2 9 


0.21 


sync 


302 


795 


408 


-0.0 


~ - +0 6 
1 


0.00 


sync 


3 


717 


924 


5.5 


-. ~+2 3 
1 Q T J 
1 ' y -3.2 


0.52 


dust 


25- - ! 


464 


216 


-0.2 


r. . + 11 

— 4 ' 


0.01 


sync 


6 


373 


786 


-0.8 


u.»_ 6 


0.00 


sync 


3 


212 


963 


1.2 


-1 -1+11 

1 - 1 -1.4 


0.28 


sync 


3 


456 


544 


-0.6 


— O-fiiiii 


0.04 


sync 


2:5 


343 


954 


2.7 


n n+1-0 
2 - 2 -1.7 


0.67 


dust 


385 


1497 


12 


1.7 


1 S+ 1 - 1 
1 - 8 -1.3 


0.52 


dust 


8 


14 


654 


-1.4 




0.11 


sync 


11 


1818 


267 


-2.5 


-l-6tl? 


0.00 


sync 


2 


1050 


183 


-0.0 


-o-i±S:i 


0.00 


sync 


4 


13 


556 


-0.7 


-0 9+ - 8 


0.00 


sync 


2 


1407 


1139 


-0.5 


-0 5+ 3 
U -°-0.3 


0.00 


sync 


3 


961 


1474 


7.3 


1 9+ 23 

-2.8 


0.52 


dust 


391 


540 


770 


7.7 


1 ' a -2.8 


0.54 


dust 


669 


55 


1086 


3.6 


1 3+ 24 
1 -°-3.2 


0.46 


sync 


268 


2115 


.511 


-0.1 


-0 2 +0 - 8 

u - z -1.0 


0.00 


sync 


6 


722 


688 


-1.2 


, , + 0.5 
± -° — 0.6 


0.00 


sync 


5 


30 


990 


-1.7 


-1 3+ 13 


0.01 


sync 


6 


577 


1578 



TABLE 5 — Continued 



ID &: coordinates: 2.0 mm data: 1.4 mm data: Spectral index &: type: Nearest source in: 

SPT ID RA DEC S/N S" 7 "" S dist S/N S raw S dist a raw a dist P(q > 1.66) Type SUMSS RASS IRAS 

[deg] [deg] [ m Jy] [ m Jy] [ m Jy] [ m Jy] [arcscc] [arcscc] [arcscc] 



SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 



J051030- 
J051040- 
J051116- 
J051121- 
J051135- 
J051140- 
J051217- 
J051254- 
J051259- 
J051355- 
J051358- 
J051400- 
J051404- 
J051406- 
J051418- 
J051424- 
J051433- 
J051439- 
J051456- 
J051459- 
J051506- 
J051517- 
J051535- 
J051550- 
J051614- 
J051618- 
J051639- 
J051659- 
J051744- 
J051811- 
J051825- 
J051841- 
J051909- 
J051939- 
J052039- 



5844.4 

5558.3 

5341. 9 C 

5920.0 

5809.5 

5846.7 

5724. a ' f 

5413.3 

5935.6 

5055.7 

5826.0 

5046.2 

5555.2 

5124.6 

5629.5 

5744.8 

5532.3 

5329.0 

5114.7 

5130.0 

5344.2 

5500.2 

5657.1 

5546.5 

5429. 6 g 

5417.8 

5920.4 

5736.9 

5851.4 

5144.1 

5614.2 

5006.1 

5609.6 

5044.8 

5329.8 



77.628 
77.670 
77.819 
77.840 
77.899 
77.919 
78.073 
78.226 
78.247 
78.483 
78.496 
78.503 
78.521 
78.527 
78.577 
78.604 
78.638 
78.666 
78.735 
78.748 
78.779 
78.824 
78.898 
78.958 
79.059 
79.077 
79.165 
79.249 
79.434 
79.547 
79.608 
79.673 
79.788 
79.914 
80.164 



-58.740 
-55.972 
-53.700 
-59.333 
-58.158 
-58.779 
-57.400 
-54.223 
-59.594 
-50.928 
-58.435 
-50.771 
-55.920 
-51.410 
-56.492 
-57.747 
-55.539 
-53.483 
-51.246 
-51.501 
-53.738 
-55.004 
-56.953 
-55.776 
-54.494 
-54.298 
-59.341 
-57.616 
-58.858 
-51.735 
-56.237 
-50.103 
-56.161 
-50.748 
-53.498 



12.43 
2.04 
5.23 
6.62 
4.72 
5.17 
4.54 
4.73 
4.69 
23.09 
5.18 
5.16 
0.92 
6.29 
7.44 
9.45 
6.61 
6.39 
4.66 
5.68 
6.05 
4.60 
16.31 
6.56 
4.67 
4.78 
1.95 
0.13 
9.13 
22.54 
16.27 
12.94 
5.11 
5.78 
3.90 



15.04 
2.56 
6.57 
8.01 
5.71 
6.26 
5.50 
5.94 
5.68 

30.27 
6.28 
6.77 

I. 16 
8.25 
9.35 

II. 44 
8.30 
8.14 
6.01 
7.32 
7.70 
5.86 
20.78 
8.36 
5.94 
6.09 
2.38 
0.15 
11.15 
29.06 
20.72 
16.68 
6.51 
7.46 
4.97 



14 77+ 1 - 35 
' ' -1.38 

1 77+ 1 - 35 
1 - ' ' -1.21 

5 91 + 133 

7 r 9 + 1.29 
'• oz -1.31 
q 71 + 3.80 

S.28il;i» 

, 7fi +2.03 
°- '"-3.36 
r. =4+1 -20 

5.46± 2 f 3 

72+ 1 05 
7.64+i;» 
8.84±i;il 
ii.oo±J;U 

7 70+ 1 ' 34 

1 ■ ,u — 1.35 

q + 1.36 
— 1.39 



7.48^ 



, ,,+4.07 
A - x — 0.89 

6 44+ 1 40 

D '^*-1.47 
+ 1.37 
-1.37 

2.02l?; 7 ° 

+ 1.52 



7.12 4 



20.51 
7.76 
2.96 
3.91 
2.55 
0.48 



- 1. 



K + 1.35 
-1.38 
K+2.60 
-2.66 
, +2.02 
-3.52 
K + l-04 
-1.40 
q + 0.72 
-0.30 

10.79±iJi 
28.80l 1 1 ; 7 4 ' 
20.46±i;|| 

d + 1.47 
8-1.44 
T+1.50 
-2.08 
K+138 
-1.44 

4.40±ij| 



16.: 

5.32^ 

6.66 4 



4.29 
4.59 
5.16 
3.89 
0.00 
0.09 
6.63 
0.58 
7.72 
8.25 
0.00 
0.52 
4.51 
1.51 
2.28 
1.92 
1.94 
1.74 
0.00 
1.04 
4.18 
0.26 
4.74 
2.18 
0.31 
0.37 
5.42 
4.60 
4.55 
7.92 
5.06 
2.24 
0.89 
1.37 
4.66 



13.74 
15.36 
17.27 
12.47 
-4.35 
0.28 
21.24 
1.94 
24.71 
28.56 
-1.31 
1.78 
15.07 
5.23 
7.61 
6.17 
6.50 
5.87 
-1.50 
3.60 
14.09 
0.87 
15.97 
7.33 
1.03 
1.26 
17.32 
14.70 
14.55 
27.44 
17.05 
7.75 
3.00 
4.75 
15.70 



12.79l 3 3 ;f 7 


-0.2 


-0 4 +0 - 7 





00 


sync 


15 


761 


891 


1.8611;?* 


4.9 


1 2+ 2 ' 7 

A - -2.8 





45 


sync 


145 


1014 


1357 


16.27±1;™ 


2.6 


„ O + 0.9 
^• 5 -0.9 





89 


dust 


203 


1025 


287 


ll.28±l;l? 


1.2 


1 1+0-9 
1 ' i -l.l 





27 


sync 


5 


742 


733 




NaN 


-i.ol?° 4 





.18 


sync 


21 


608 


881 


2 64+ 1 69 
z -°*-1.16 


-8.5 







06 


sync 


13 


439 


459 


19.l8t|;™ 


3.7 


Q c+0.8 
J -°-0.8 





98 


dust 


118 


765 


4 


2.36l?; 7 | 


-3.0 


-0.8l 2 j 





16 


sync 


18 


259 


414 


22.73±|:il 


4.0 


3 9+ 06 


1 


00 


dust 


204 


540 


162 


27.78±1;|I 


-0.2 


-0 2+°- 4 





00 


sync 


2 


578 


1047 


2 39+ 1 ' 44 

iiJa - 1 .25 


NaN 


-1 7+ 17 

1 - ' -1.0 





06 


sync 


22 


262 


673 


„ , c+2.13 
°- io -1.31 


-3.6 


-1 3+ 17 

— 1.2 





06 


sync 


24 


8 


491 


1 06 +lsl 


7.0 


1 7+2.3 





51 


dust 


282 


540 


1119 


S.ool 2 ; 8 .* 


-1.2 







01 


sync 


5 


1148 


197 


K fic; + 3 14 

6.65_ 2 47 


-0.6 


-o.sll; 2 





01 


sync 


1 


728 


390 


6 21+ 2 ™ 


-1.7 


-i-6lo:°9 





00 


sync 


4 


362 


1007 


5 69 +2 96 

d - u -2.16 


-0.7 


-o.slij 





01 


sync 


10 


449 


566 


e: on + 2.87 
-1.98 


-0.9 


-0-9111 





01 


sync 


23 


1436 


941 


i.36ti:S 


NaN 







19 


sync 


28 


1797 


571 


4.05l 2 ;f 5 


-1.9 


-1 2+ 14 

- L-Z — 1.2 





02 


sync 


5 


990 


668 


12.93l|; 7 § 


1.6 


1 6 +0.9 
1 ' D — 1.0 





-IS 


sync 


8 


971 


12 


L - ' a -1.44 


-5.2 


-0.71 2 ? 





20 


sync 


275 


1168 


347 


ls.oeli;^ 


-0.7 


-o.sir, 





00 


sync 


15 


1865 


959 


6 21+ 3 ' 04 

u --- L -2.44 


-0.4 


-0-611; 2 





01 


sync 


4 


773 


1328 


1 QC + 1 " 
1 - a °-1.60 


-4.8 


-O.sl 2 ; 8 





.18 


sync 


104 


193 


783 


„ ,0+2.01 
.77 


-4.3 


-o.9l 2 ; 4 4 





15 


sync 


12 


505 


943 


9 86 +5 ' 75 


5.4 


3 6+ 10 





75 


dust 


189 


266 


897 


90+ 1 38 
u - au -0.52 


12.4 


2.2l 20 7 





58 


dust 


253 


257 


763 


13.54l 3 3 ;f 5 


0.7 


6+°- 7 





06 


sync 


5 


251 


1436 


26.6711;?° 


-0.2 


-0 2+°' 4 





00 


sync 


3 


652 


665 


I6.i4l 3 3 ; 7 7 


-0.5 


-0 6+°' 6 





00 


sync 


13 


95 


650 


o 1 £ + 2-88 
°- i0 -1.91 


-2.1 


-1.9l«; 8 7 





00 


sync 


14 


634 


423 


o OJ-+2.34 
0,0 ° — 1.43 


-2.1 


-1 0+ 17 





06 


sync 


4 


497 


714 


4 54+ 2 ' 71 

^■ J ^-1.71 


-1.2 


-i-olll 





112 


sync 


9 


7 


532 


S qo+6.77 
S - JS -6.30 


3.1 


i-6l^; 8 5 





50 


sync 


297 


1078 


865 



TABLE 5 — Continued 



ID &: coordinates: 2.0 mm data: 1.4 mm data: Spectral index &: type: Nearest source in: 

SPT ID RA DEC S/N S" 7 "" S dist S/N S raw S dist a raw a dist P(q > 1.66) Type SUMSS RASS IRAS 

[deg] [deg] [ m Jy] [ m Jy] [ m Jy] [ m Jy] [arcscc] [arcscc] [arcscc] 



SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 



J052044- 
J052124- 
J052137- 
J052206- 
J052215- 
J052309- 
J052318- 
J052400- 
J052405- 
J052440- 
J052500- 
J052558- 
J052731- 
J052743- 
J052759- 
J052832- 
J052834- 
J052841- 
J052846- 
J052850- 
J052852- 
J052903- 
J052904- 
J053058- 
J053107- 
J053117- 
J053147- 
J053205- 
J053208- 
J053244- 
J053248- 
J053250- 
J053310- 
J053312- 
J053319- 



5508.4 

5412.5 

5456.7 

5016.3 

5607.4 

5822.9 

5308.5 

5133.9 

5602.2 

5658.8 

5135.2 

5725.1 

5546.0 

5426.3 

5232.2 

5433.8 

5820. 3 a 

5726.0 

5919.9 

5300. 3 h 

5457.6 

5436.6 

5538.8 

5951.3 

5543.5 

5504.4 

5414.4 

5434.6 

5310.5 

5812.1 

5721.6 

5047.1 

5453.3 

5238. 5 a 

5022.0 



80.187 
80.353 
80.407 
80.528 
80.565 
80.789 
80.826 
81.003 
81.021 
81.168 
81.250 
81.495 
81.883 
81.931 
81.998 
82.137 
82.145 
82.174 
82.192 
82.212 
82.220 
82.263 
82.271 
82.743 
82.780 
82.822 
82.946 
83.024 
83.035 
83.186 
83.201 
83.212 
83.293 
83.302 
83.332 



-55.141 
-54.210 
-54.946 
-50.272 
-56.124 
-58.383 
-53.143 
-51.566 
-56.037 
-56.981 
-51.588 
-57.419 
-55.767 
-54.438 
-52.537 
-54.564 
-58.339 
-57.433 
-59.333 
-53.006 
-54.960 
-54.611 
-55.647 
-59.856 
-55.725 
-55.074 
-54.241 
-54.577 
-53.176 
-58.202 
-57.360 
-50.786 
-54.889 
-52.642 
-50.367 



16.69 
10.84 
8.67 
4.59 
5.24 

I. 96 
9.18 
7.74 

II. 32 
106.29 

8.09 
54.06 
1.84 
16.52 
2.49 
7.53 
19.89 
1.93 
12.49 
1.13 
6.87 
7.34 
6.44 
12.66 
8.76 
28.10 
6.66 
1.09 
41.98 
6.02 
4.91 
10.53 
4.94 
1.65 
5.33 



21.26 
13.81 
11.05 
5.92 
6.68 
2.40 
11.69 
9.97 
14.41 
135.38 
10.42 
66.02 
2.34 
21.04 
3.17 
9.59 
24.29 
2.35 
15.25 

I. 44 
8.75 
9.34 
8.21 
15.46 

II. 16 
35.80 
8.49 
1.39 

53.48 
7.35 
6.00 
13.58 
6.29 
2.11 
6.87 



2i.ooti;fi 
13.45±i;S 
io.65±i;|l 

. ,,, + 1.43 

I 67+ 1 * 30 

1DI -1.13 

II i n + 1 - 38 
11 - lz -1.36 
a rn + 1-39 

14.07±i;^ 
134.67^;^ 

Q Q4+ 1 ' 37 

65 " 67 -2.97 
A ,<l ' -1.02 

20.77± 1 1 ;f 1 
2.48t 1 1 ;f 7 

n itn + 136 

a.io_ 1 36 



24.07t;^ 
1.65 



+ 1.30 
' — 1.11 

14.99±—® 

1 11+ 1 - 41 
-'-•- L1 -o.so 

Q+1.36 
y -1.36 

9 17+ 131 

K - x 1 -1.31 

d + 1.36 
3 -1.38 

35.54±-SI 

•3+137 
-1.37 
, 19 + 1.43 

1 -- Li -o.si 
53.19t^; 5 5 i 



8.191 



10.68] 



7.73 4 



+ 1.31 
-1.34 

c + 1-43 
-1.87 

13.35±1;^ 

1 +1.41 
-1.43 

•3+137 
-1.21 



6.727 
4.96j 
13.35 
5.5V 
1.73 4 



r nK + 1.42 



6.54 
2.74 
3.65 
1.44 
1.50 
4.61 
1.25 
3.74 
3.29 

32.01 
2.21 

16.50 
4.51 
6.83 
4.65 
4.80 
6.63 
4.62 
5.47 
5.09 
2.41 

10.96 
0.00 
5.67 
2.02 
7.84 
0.79 
5.14 

14.79 
2.23 
1.97 

12.20 
4.12 
5.04 
1.47 



22.02 
9.23 
12.28 
4.98 
5.07 
14.74 
4.21 
12.97 
11.07 
107.78 
7.64 
52.71 
15.18 
23.01 
15.67 
16.18 
21.17 
14.76 
17.47 
17.14 
8.11 
36.89 
-0.02 
18.11 
6.78 
26.40 
2.67 
17.30 
49.78 
7.14 
6.31 
42.26 
13.88 
16.97 
5.09 



21.20 
8.56 



3.96 
3.92 
+ 3.26 



n.02t 3 3 ; 4 5 8 5 



-, + 2.85 
-2.21 
+ 2.69 
1.80 
+6.85 
1.16 
+ 2.32 
1.36 
+ 3.62 
3.72 



3.12 
4.28 

I. 87 
5.55 

II. 66 
10.12 

106.49 

6.89 
51.93i= 

1.58 
22.24 

2.75 
15.18 
20.42 

1 



+ 3.49 
-3.14 
+ 10.15 
10.15 
K3.22 
-2.42 



+ 4.90 
0.97 
+ 3.97 
4.01 
+9.52 
1.76 
+ 3.67 
3.79 
+ 3.73 
3.78 



1+6.72 
— 1.13 

'"■ w -a .64 

+8.06 
1.06 
+ 3.12 
-2.67 
+4.77 
4.77 
+ 1.66 
-0.94 

I 7 q-7+3.63 

II - zl -3.63 
,+2.97 

2.05 
+4.09 
4.16 
+ 2.16 
1.23 
+8.40 



1.59 
6.85 

35.38 
3.56 



6.52 
25.67 
4.18 
1.61 



48.98± B 6 ; 6 5 ? 



5.81 
4.56 
40.77 
12.61 
2.34 
4.39 



+ 2.77 
2.38 



+ 2.71 
2.22 
+ 5.23 
5.23 
3.59 
3.75 
+ 10.12 

1.62 
+ 2.74 
1.82 



0.1 


o o +0 - 5 


0.00 


sync 


3 


1012 


1653 


-1.1 


1 n + 0.9 
1 ■ z - 1 . 


0.00 


sync 


9 


1109 


128 


0.3 


o-ili.! 


0.01 


sync 


8 


514 


829 


-0.5 




0.20 


sync 


395 


784 


592 


-0.8 


-0.7±i;l 


0.06 


sync 


13 


676 


757 


4.9 


1 4+ 2 ' 6 

— 2.9 


0.47 


sync 


21!) 


995 


376 


-2.8 


_1 Q+io 

— 0.8 


0.00 


sync 


1 


661 


1064 


0.7 


0.6l 1 ^ 


0.08 


sync 


7 


503 


1247 


-0.7 


- .9±r 


0.00 


sync 


(i 


72.S 


1083 


-0.6 


-0.6l° o j 


0.00 


sync 


2 


1078 


798 


-0.8 


-1 0+ 11 

l.U_ 1 2 


0.00 


sync 


5 


79 


1217 


-0.6 


-0.6l° j 


0.00 


sync 


2 


789 


347 


5.1 


i.it 2 2 s 7 


0.44 


sync 


128 


1413 


1607 


0.2 


2+°° 

u - -0.6 


0.00 


sync 


5 


406 


1340 


4.4 


1 o + 2.6 
A '°— 3 


0.47 


sync 


276 


1773 


1001 


1.4 




0.35 


sync 


13 


269 


1582 


-0.4 


-0.5t° 5 5 


0.00 


sync 


23 


409 


131 


5.0 


l-4± 2 2 ; 6 9 


0.47 


sync 


190 


1438 


442 


0.4 


0.3l° 6 7 


0.00 


sync 


(i 


7 


708 


0.7 


9 O + 2.0 
Z.D_3 2 


0.58 


dust 


45 


850 


Oil 


-0.2 


-0.411:2 


0.01 


sync 


30 


1278 


1317 


3.7 

NaN 


q 7+05 
°' ' -0.5 

-2 0+ 11 

• £ - u — .7 


1.00 
0.00 


dust 
sync 


121 
4 


543 
1165 


1757 
1122 


0.4 


4 +0 - 6 


0.01 


sync 


(i 


1326 


910 


-1.4 


-1 4+ 11 


0.00 


sync 


7 


1320 


272 


-0.8 


-0-9lo1 


0.00 


sync 


4 


33 


752 


-3.1 




0.00 


sync 


8 


81)2 


412 


6.9 




0.59 


dust 


493 


263 


343 


-0.2 


-0-2l j 


0.00 


sync 


2 


73 


148 


-0.1 


-0.31- 


0.03 


sync 


9 


604 


1034 


0.1 


1 + 1 6 


0.13 


sync 


330 


375 


(ill 


3.1 
2.2 


3 n+o.4 

2.3±i : ; 


1.00 
0.71 


dust 
dust 


424 

8 


1002 
736 


441 

300 


5.7 


2 4+ 20 

^•^-3.4 


0.58 


dust 


12 


744 


9 


-0.8 


-0 7+ 15 

u - ' -1.5 


0.04 


sync 


8 


597 


783 



TABLE 5 — Continued 



ID &: coordinates: 2.0 mm data: 1.4 mm data: Spectral index &: type: Nearest source in: 

SPT ID RA DEC S/N S" 7 "" S dist S/N S raw S dist a raw a dist P(q > 1.66) Type SUMSS RASS IRAS 

[deg] [deg] [ m Jy] [ m Jy] [ m Jy] [ m Jy] [arcscc] [arcscc] [arcscc] 



SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 



J053324- 
J053327- 
J053345- 
J053412- 
J053426- 
J053458- 
J053527- 
J053631- 
J053644- 
J053646- 
J053726- 
J053748- 
J053816- 
J053819- 
J053823- 
J053834- 
J053849- 
J053857- 
J053909- 
J053942- 
J054020- 
J054025- 
J054030- 
J054045- 
J054120- 
J054123- 
J054135- 
J054223- 
J054254- 
J054343- 
J054357- 
J054407- 
J054431- 
J054434- 
J054510- 



5549.6 

5050.8 

5818. 1" 1 

5924.3 

5113.2 

5439.0 

5013.6 

5258.8 

5310.5 

5714.1 

5434. 4> 

5718.4 

5030.8 

5227.7 

5625.6 

5911.0 

5955.5 

5712.2 

5511.0" 

5612. 8 a 

5356.0 

5303.7 

5356.5 

5418.3 

5738.3 

5752.8 

5016.9 

5142.9 

5122.0 

5813.3 

5532.1 

5444.2 

5826.2 

5402.9 

5635.1 



83.352 
83.366 
83.441 
83.553 
83.610 
83.744 
83.864 
84.130 
84.185 
84.194 
84.360 
84.453 
84.569 
84.580 
84.598 
84.643 
84.706 
84.738 
84.789 
84.927 
85.085 
85.105 
85.126 
85.191 
85.337 
85.348 
85.400 
85.598 
85.728 
85.931 
85.991 
86.032 
86.131 
86.144 
86.294 



-55.827 
-50.847 
-58.302 
-59.406 
-51.221 
-54.651 
-50.227 
-52.981 
-53.176 
-57.236 
-54.574 
-57.308 
-50.514 
-52.462 
-56.427 
-59.184 
-59.926 
-57.204 
-55.183 
-56.214 
-53.935 
-53.063 
-53.942 
-54.306 
-57.640 
-57.882 
-50.282 
-51.716 
-51.368 
-58.223 
-55.536 
-54.737 
-58.438 
-54.049 
-56.586 



55.02 
7.82 
7.24 
15.34 
10.88 
20.37 
0.62 
1.18 
1.15 
12.01 
6.36 
12.53 
6.81 
8.32 
5.22 
6.31 
4.07 
4.72 
17.25 
99.72 
4.59 
23.75 
15.42 
259.54 
2.24 
4.27 
4.98 
70.43 
2.48 
4.67 
14.32 
0.00 
5.88 
0.21 
0.51 



70.08 
10.09 
8.84 
18.73 
14.02 
25.95 
0.79 
1.50 
1.46 
15.30 
8.11 
15.96 
8.78 
10.60 
6.64 
7.70 
4.97 
6.01 
21.97 
127.03 
5.84 
30.26 
19.64 
330.59 
2.73 
5.21 
6.41 
90.80 
3.20 
5.70 
18.24 
-0.41 
7.18 
0.27 
0.64 



69.71±3 }s 
9.48+.^ 

„ ,3+1.32 
0-50-1 .30 

18.46±1; 44 
l3.6l±J;g 

26.71+i;Sjj 
0.65ii;S5 

n Q4+ 1 - 28 

u ' y4: -0.66 

0.90+. ; 24 
I4.99±i;tl 
7.36±1; 3 | 

15.66±i^ 
8.54+i : |l 
10.19t\l 5 7 
5.78tl;S 
r.i6+i;U 

i.47±\f a 
4 57+ 1 66 

*•»>' -3.40 

21.69±i;j£ 
126.35+. 5 5 ; 4 J 

2.071?;- 
30.04tl;2 
19.3511;" 

+ 13.64 
-13.02 
+ 1.18 
1.39 
+ 1.26 
1.29 
+ 1.49 



-1. 



328 
2.35 
4.74 
5.36 
90.32+. 3 ;;) 2 , 

,1 + 1-41 
-1.66 

17.9715;*? 

, +0.64 
-0.27 

R+1.31 

-1.36 
1 +0.75 
-0.35 
n +0.88 
-0.41 



4.743 



0.41: 

6.46^ 

0.51; 

0.60j 



17.54 
1.35 
3.66 
3.96 
2.15 
7.12 
4.87 
4.72 
4.67 
3.92 
1.07 
5.32 
9.04 
3.68 
2.10 
3.59 
4.74 
1.67 
3.69 

29.04 
0.33 
6.57 
3.86 

78.29 
4.93 
4.78 
1.84 

20.28 
4.50 
3.02 
4.56 
4.83 
1.53 
4.57 
4.62 



59.06 

4.69 
11.71 
12.64 

7.44 
23.96 
16.88 
15.89 
15.71 
13.20 

3.61 
17.91 
31.31 
12.40 

7.08 
11.46 
15.14 

5.61 
12.42 
97.77 

1.11 
22.13 
12.98 
263.58 
15.74 
15.28 

6.38 
70.27 
15.61 

9.63 
15.34 
16.27 

4.87 
15.40 
15.57 



58.28 
5.32 



+ 6.18 
6.31 
+ 2.59 



10.48t 3 3 ;^ 



11.85 
7.46 



+ 3.37 
3.13 
+ 3.00 
2.02 



+ 2.09 
0.66 
+ 3.28 
0.77 



23.18+T t 



1.06 
1.27 

I. 22 
12.1 

4.32 
17.01 
29.68 

II. 15 
5.41 

10.17 
8.96 
3.77 

11. 



+ 2.77 
-0.77 
+ 3.61 
-3.53 
t-2.40 
-1.40 
+ 3.73 
3.83 
+4.52 
-4.58 
3.49 
3.58 
+ 2.87 
2.44 
+ 3.33 
3.46 
+ 5.91 
6.70 
+ 2.78 
-2.12 
+ 3.38 
-2.88 



96.58+lf 9 



1.79 



+ 2.03 
1.47 

21.36+.3j£ 

i n -i<7+3.53 
' -3.20 

260.65l 2 2 4 ;^ 

,+9.35 



3.73 
9.75 
4.62 
69.38 
2.24 
7.83 



14.38t 3 ;^ 



2.71 
+ 5.50 
7.25 
+ 2.78 
2.19 
+ 7.15 
7.15 
+8.66 
1.35 
+ 3.32 
3.56 



0.81 
4.49 

0.88 
0.95 



+ 1.25 
0.46 
+ 2.60 
1.70 
+ 1.33 
0.52 
+ 1.58 
0.55 



-0.5 


-0 5 +0 - 3 

U - J — 0.3 


0.00 


sync 


5 


1224 


975 


-2.1 


-1 6+ 11 


0.00 


sync 


12 


793 


392 


0.8 




0.10 


sync 


36 


4 


514 


-1.1 


-1 9 + 0.7 
x - z — 0.8 


0.00 


sync 


5 


49 


1561 


-1.7 


-1 6 +0 - 9 


0.00 


sync 


1 


473 


742 


-0.2 


-0-3-0.1 


0.00 


sync 


3 


570 


907 


8.3 


2.0t 2 2 2 s 


0.55 


dust 


30 


645 


419 


6.4 




0.51 


dust 


167 


865 


252 


6.5 


l-7± 2 2-.8 


0.51 


dust 


483 


951 


690 


-0.4 


-0 6 +0 - 8 
u -°— 1.0 


0.00 


sync 


3 


353 


826 


-2.2 


-1 5+ 13 
± -° — 1.0 


0.00 


sync 


83 


1366 


192 


0.3 


2+ - 6 

u - -0.7 


0.00 


sync 





16 


884 


.3.5 


3 4 +0 - 6 

° -0.6 


1.00 


dust 


190 


559 


454 


0.4 


0.3+.™ 


0.03 


sync 


10 


3 


409 


0.2 


o.oli 1 


0.09 


sync 


14 


648 


669 


1.1 


i.oii-3 


0.23 


sync 


14 


414 


275 


3.0 


1.8ii;« 


0.51 


dust 


609 


142 


2:52 


-0.2 


-o.oi; : g 


0.16 


sync 


14 


490 


1361 


-1.6 


-1 6+ - 7 

-■-■ D -0.8 


0.00 


sync 


8 


843 


832 


-0.7 


-0 7 +0 - 2 

u - ' -0.3 


0.00 


sync 


5 


689 


1244 


-4.5 


-o-6l?:? 


0.21 


sync 


19 


759 


1145 


-0.9 


-0.9t°i 


0.00 


sync 


2 


503 


643 


-1.1 


-1 3+ 07 
J-- — 8 


0.00 


sync 


5 


810 


1226 


-0.6 


-o-e±S:3 


0.00 


sync 


5 


56 


191 


4.8 




0.57 


dust 


294 


657 


461 


2.<> 


1 R+ 1 - 5 
1-8-3.4 


0.52 


dust 


73 


798 


1293 


-0.0 


-o.it\- 6 7 


0.10 


sync 


91 


1277 


1082 


-0.7 


-o.7+_° : i 


0.00 


sync 


4 


9 


914 


4.3 


0.8l 3 ;? 


0.41 


sync 


408 


187 


853 


1.4 


i -°-1.6 


0.43 


sync 


1!) 


469 


375 


-0.5 


-0 6+ ' 6 


0.00 


sync 


10 


6 


155 


NaN 




0.61 


dust 


405 


531 


800 


-1.1 


-o.9li:S 


0.02 


sync 


10 


416 


553 


11.0 




0.57 


dust 


622 


694 


851 


8.7 


z ' u -2.8 


0.54 


dust 


428 


510 


1406 



TABLE 5 — Continued 



ID &£ coordinates: 2.0 mm data: 1.4 mm data: Spectral index &: type: Nearest source in: 

SPT ID RA DEC S/N 5 raw 5 dist S/N 5 raw 5 diat a raw a dist P(a > 1.66) Type SUMSS RASS IRAS 

[deg] [deg] [niJy] [niJy] [mJy] [ m Jy] [arcscc] [arcscc] [arcsec] 



SPT- 


■S 


J054545- 


-5138 


.0 


86.439 


-51 


.635 


SPT- 


■S 


J054554- 


■5815 


.0 


86.478 


-58 


251 


SPT- 


■S 


J054644- 


-5858 


,3 


86.684 


-58 


.973 


SPT- 


■S 


J054716- 


-5104 


2 k 


86.818 


-51 


.070 


SPT- 


■S 


J054733- 


-5335 


.6 


86.888 


-53 


.594 


SPT- 


■S 


J054804- 


-5955 


.0 


87.020 


-59 


917 


SPT- 


■S 


J054830- 


-5218 


.5 


87.126 


-52 


.308 


SPT- 


■S 


J054846- 


-5748 


,8 


87.194 


-57 


.814 


SPT- 


■S 


J054901- 


-5752 


.4 


87.258 


-57 


.874 


SPT- 


■S 


J054903- 


-5741 


3 


87.265 


-57 


.689 


SPT- 


■S 


J054912- 


-5026 


.6 


87.304 


-50 


.444 


SPT- 


■S 


J054941- 


-5645 


.0 


87.422 


-56 


.751 


SPT- 


•s 


J054943- 


-5246 


.4 


87.432 


-52 


.774 


SPT- 


■s 


J054951- 


-5047 


.8 


87.465 


-50 


.798 


SPT- 


■S 


J054953- 


-5358 


.4 


87.472 


-53 


.975 


SPT- 


■S 


J055002- 


-5356 


.6 


87.509 


-53 


.943 


SPT- 


■S 


J055009- 


-5732 


3 


87.540 


-57 


.540 


SPT- 


■s 


J055046- 


-5304 


.9 


87.694 


-53 


.082 


SPT- 


■S 


J055116- 


-5334 


.4 1 


87.817 


-53 


574 


SPT- 


•s 


J055117- 


-5007 


.4 


87.823 


-50 


124 


SPT- 


■s 


J055119- 


-5545 


.5 


87.830 


-55 


.760 


SPT- 


■S 


J055133- 


-5655 


,3 


87.890 


-56 


.923 


SPT- 


■s 


J055135- 


-5902 


.7 


87.897 


-59 


.046 


SPT- 


■S 


J055138- 


-5058 


.0 


87.912 


-50 


.968 


SPT- 


■s 


J055152- 


-5526 


,5 a 


87.971 


-55 


.442 


SPT- 


■s 


J055201- 


-5951 


.4 


88.005 


-59 


.858 


SPT- 


■S 


J055232- 


-5349 


4 a 


88.135 


-53 


.823 


SPT- 


•s 


J055233- 


-5242 


.7 


88.139 


-52 


712 


SPT- 


■S 


J055241- 


-5238 


.5 


88.172 


-52 


.643 


SPT- 


•s 


J055302- 


-5548 


.7 


88.259 


-55 


812 


SPT- 


■S 


J055320- 


-5007 


.3 


88.335 


-50 


.122 


SPT- 


■s 


J055331- 


-5131 


.7 


88.380 


-51 


.529 


SPT- 


•s 


J055359- 


-5051 


.8 


88.498 


-50 


.864 


SPT- 


■S 


J055421- 


-5018 


.6 


88.588 


-50 


311 


SPT- 


■s 


J055422- 


-5304 


2 


88.596 


-53 


.071 



7.85 
4.58 
1.59 
7.40 
0.01 
6.15 

19.55 
0.79 
5.89 
7.92 
9.19 
1.03 

163.42 
5.58 
4.59 
2.69 
523.43 

17.12 
4.83 
2.29 
4.59 
4.60 
6.35 
2.91 

27.63 
0.49 
9.99 

I. 61 
0.81 
3.64 
2.69 
0.96 
4.80 

II. 08 
2.20 



10.12 
5.59 
1.94 
9.75 
0.02 
7.67 

25.67 
0.99 
7.34 
9.87 

12.11 
1.35 
214.59 
7.35 
6.03 
3.54 
652.66 

22.48 
6.34 
3.01 
6.03 
6.03 
7.92 
3.84 

36.28 
0.61 

13.12 
2.12 
1.06 
4.79 
3.54 
1.26 
6.33 

14.60 
2.88 



9.33+i;» 

2.16^.2 

-19+1.13 
'■ A ^-1.43 

Q r;o + l-38 
i.3g 

o.«tg;g 
6.89+i;|S 

25.411}- 
0.68+S;"? 

a C t +1.34 
6 - 61 -1.3S 

9-^—1.33 

ii.68+i;£ 

213.45tt;t| 

6 - 50 -1.50 

5 17+ 1 - 44 
-1.47 

, „„+1.02 
3.»8_ 091 

649 19+ 26 - 92 
22.24+};^ 

2.26il;S 

J - ' ' -3.37 

7 69+ 1 29 

1 29 



r r.o + 0.94 

o.uo_ 76 
36.02tjg 7 
<\7+ - 87 

V DI -0.39 



12.73±i;J? 



1.27 
0.68 
4.21 
2.95 
0.82 
4.91 
14.26 
2.34 



+ 1.41 
0.91 
+ 1.04 
0.46 
+ 1.38 
1.32 
+ 1.33 
1.51 
+ 1.24 
0.59 
+ 1.67 
■3.44 
+ 1.46 
1.46 
+ 1.33 
1.52 



0.10 
0.43 
5.41 
7.86 
4.96 
0.93 
8.19 
4.52 
1.61 
3.67 
2.51 
4.52 

51.06 
1.47 
4.24 
6.15 

173.64 
5.37 
6.87 
4.65 
1.26 
1.11 
8.73 
9.05 
9.33 
4.56 
2.78 
4.57 
4.51 
4.78 
4.76 
4.78 
1.44 
4.67 
4.83 



0.34 

1.36 
17.28 
27.22 
16.96 

3.06 
27.97 
14.80 

5.28 
12.01 

8.70 
15.45 
174.45 

5.08 
14.48 
21.02 
568.60 
18.36 
23.48 
16.11 

4.31 

3.80 
28.59 
31.36 
31.86 
14.93 

9.48 
15.62 
15.42 
16.33 
16.48 
16.57 

5.00 
16.19 
16.52 



+ 1-69 
0.97 
+ 1.97 



4.30 
1.80 
6.21 
25.49 
0.91 
3.99 
27.24 
1.04 
4.73 
10.77 
7.91 
1.12 
172.41 
4.61 
13.15 
17.28 
562.31 
17.47 
21.46 
2.34 
2.88 
2.69 
26.99 
26.73 
31.12 
0.94 
8.65 
1.45 
1.03 
9.75 
4.25 
1.16 
3.70 



2.93t 2 ° f 



1.45 
+ 8.26 
-5.25 
+4.39 
-4.39 
+ 1.59 
-0.54 
2.26 
-1.28 
+4.25 
4.30 
+ 1.72 
-0.64 
+ 2.71 
1.82 
+ 3.40 
3.49 
+3.34 
2.63 
+ 1.97 
0.70 
+ 16.09 
15.69 
+ 2.77 
1.80 
+ 3.66 
3.85 
+4.02 
4.49 
+ 51.40 
50.17 
+3.84 
3.84 
+4.09 
-4.23 
+9.52 
1.48 
+ 2.76 
2.15 
+ 2.65 
2.06 
+4.23 
4.29 
+4.23 
4.23 
+4.49 
4.49 
+ 1.54 
0.56 
+ 3.40 
2.80 
+4.42 
0.91 
+ 1.68 
0.61 
+6.21 
7.65 
+ 10.04 

3.03 
+ 2.76 
0.71 
+ 2.80 
1.97 



-9.3 


-2 2+ 10 

-^-^ — 0.6 


0.00 


sync 


10 


38 


819 


-3.8 


-0.6±i 7 


0.20 


sync 


272 


508 


420 


6.0 
2.8 


9 7+0.6 
z - 1 -0.6 


0.71 
0.94 


dust 
dust 


301 
350 


509 
1255 


1209 
11 


19.0 


2.3± 2 3 .° 


0.59 


dust 


000 


72.3 


.322 


-2.5 


-i-5±i; 3 o 


0.00 


sync 


8 


403 


433 


0.2 


2+ 04 

u - -0.5 


0.00 


sync 


3 


457 


455 


7.4 


1 a + 2.3 
i -°-2 8 


0.52 


dust 


214 


916 


1032 


-0.9 


-o-9±i:S 


0.02 


sync 


5 


72.S 


020 


0.5 


0At°i 


0.05 


sync 


2 


700 


839 


-0.9 


-1.111;? 


0.00 


sync 


10 


2.S 


591 


6.6 


1 6+ 2 ' 4 


0.49 


sync 


221 


235 


1029 


-0.6 


-0 6 +0 ' 2 

u- -0.2 


0.00 


sync 


4 


22 


369 


-1.0 


-0.91H 


0.03 


sync 


15 


530 


070 


2.4 


2.611 1 


0.79 


dust 


ISO 


1149 


1087 


4.9 
-0.4 


4-i±S:S 

-0 4+°- 2 

u- -0.2 


0.96 
0.00 


dust 

sync 


101 

2 


1069 
(i 


951 
624 


-0.6 


_n 7+0-6 

u- ' -0.7 


0.00 


sync 


3 


211 


879 


3.6 


3 4+ ' 8 

° • ^ - . 8 


0.98 


dust 


14 


732 


4 


4.6 


1.3l^ 


0.47 


sync 


207 


636 


480 


-0.9 


-0.2111 


0.19 


sync 


8 


911 


1645 


-1.3 


-0-3i?; 5 8 


0.20 


sync 


14 


457 


190 


3.5 
5.7 
-0.4 


3 4+ 06 

J " -0.6 

4 6+ - 3 

-o.4ir 4 


1.00 
1.00 
0.00 


dust 
dust 
sync 


7 

260 
4 


10.3 
1281 
550 


4 
70 
733 


8.7 


2-0l 2 2 i 


0.55 


dust 


255 


226 


124 


-0.9 


-1 0+ 09 
J - u — i.i 


0.00 


sync 


5 


526 


071 


5.4 


1.3l 2 2 ; 6 8 


0.46 


sync 


464 


647 


1162 


7.3 


1.7± 2 2- 3 8 


0.51 


dust 


487 


756 


1155 


3.3 


2 1 + 1 6 


0.56 


dust 


52!) 


667 


1820 


4.2 


1 q+ 2 - 2 


0.52 


dust 


543 


937 


1002 


7.0 


1 9+ 23 


0.52 


dust 


263 


341 


.351 


-0.6 


-0.311;? 


0.13 


sync 


4 


98 


701 


0.3 


2+°' 7 


0.00 


sync 


4 


658 


217 


4.7 


2 0"*" 2 ' 2 
z.u_ 3 4 


0.53 


dust 


327 


274 


1578 



TABLE 5 — Continued 



ID &: coordinates: 2.0 mm data: 1.4 mm data: Spectral index &: type: Nearest source in: 

SPT ID RA DEC S/N S" 7 "" S dist S/N S raw S dist a raw a dist P(q > 1.66) Type SUMSS RASS IRAS 

[deg] [deg] [mJy] [mJy] [ m Jy] [ m Jy] [arcsec] [arcscc] [arcscc] 



SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 
SPT-S 



J055455- 
J055512- 
J055528- 
J055545- 
J055555- 
J055701- 
J055732- 
J055734- 
J055811- 
J055830- 
J055833- 
J055854- 
J055909- 
J055923- 
J055947- 



5541.3 

5044.9 

5154.3 

5605.1 

5521.8 

5902.6 

5359. 3 a 

5801.0 

5029.8 

5326.4 

5835.5 

5333.8 

5128.1 

5900.6 

5026.8 



88.731 
88.802 
88.868 
88.939 
88.979 
89.258 
89.384 
89.394 
89.549 
89.627 
89.638 
89.728 
89.788 
89.846 
89.947 



-55.690 
-50.750 
-51.906 
-56.085 
-55.364 
-59.044 
-53.988 
-58.018 
-50.497 
-53.440 
-58.593 
-53.564 
-51.469 
-59.010 
-50.447 



5.91 
6.98 
5.00 
5.00 
4.65 
4.97 
5.44 
6.37 

50.99 
8.83 
6.43 
5.01 
1.25 
5.70 

37.59 



7.77 
9.20 
6.59 
6.56 
6.10 
6.19 
7.15 
7.94 
67.18 
11.59 
8.02 
6.58 
1.64 
7.10 
49.52 



7 04+ 1 41 
4 89+ 1 74 

'■° 3 -4 .14 
5-19^2.88 

6 43+ 1 41 

°-*°-1.49 
- ,9+1.34 
' ""-1.34 

66.83j:3'g2 

n.oi - l 4 ', 

9 97+1 ■ 18 
-1.66 

49.26 4 ^' 41 1 



1.85 
1.85 
0.22 
0.87 
0.00 
1.61 
2.83 
2.04 
16.10 
1.35 
2.63 
3.91 
5.68 
4.26 
12.09 



6.30 
6.40 
0.76 
2.99 
-0.34 
5.26 
9.67 
6.69 
55.77 
4.60 
8.63 
13.35 
19.68 
13.94 
41.90 



5.38 
5.86 
2.70 
3.25 
1.57 
4.02 
8.00 
5.71 



+ 2.93 
2.13 
+3.04 
-2.09 
+ 1.94 
1.62 
+ 2.39 
1.54 
+ 1.86 
1.30 
+ 2.66 
1.89 
+ 3.35 
3.37 
+ 2.90 
2.25 



54 98+ 6 05 



5.66 
7.16 4 



+ 2.46 
1.46 



19 ni +3 63 
±Z.U1_ 3 gl 



8.54 
12.82 
41.1 



+ 7.82 
7.64 

-3.50 
-3.66 
+ 5.11 
-5.11 



-0.6 
-1.0 
-5.9 
-2.1 
NaN 
-0.4 
0.8 
-0.5 
-0.5 
-2.5 
0.2 
1.0 
0.8 
1.8 
-0.5 



+ 1.3 
-1.4 
-1+1.2 
-1.2 
+ 2.0 
-1.2 
-.+ 1.8 
J — 1.4 
j + 3.0 
-1.6 
+ 1.7 
-1.7 
7+12 
-1.5 
+ 1.3 
-1.4 
- + 0.3 
J -0.3 
j+1.0 
-0.8 
+ 1.2 
-1.4 

1+1-1 
-1.2 
= +0.9 
-3.5 
j+0.9 
' — 1.0 



-0.7 
-1.0 
-1.3 
-1.0 
-0.8 
-0.3 

0.7 
-0.6 
-0.5 
-1.8 
-0.0 

2.0 

3 

1.9 

-0 5+ 3 



0.02 
0.00 
0.09 
0.08 
0.20 
0.10 
0.18 
0.01 
0.00 
0.00 
0.04 
0.62 
0.75 
0.58 
0.00 



sync 
sync 
sync 
sync 
sync 
sync 
sync 
sync 
sync 
sync 
sync 
dust 
dust 
dust 
sync 



6 
11 

9 

52 
17 

5 

10 
15 

1 



9 
11 
750 

(i 

5 



90 
803 
790 
398 
540 
295 
192 
1347 
925 
383 
401 
455 
1094 
514 

12 



1353 
603 
717 
501 
810 
546 
202 
002 
137 
365 

1445 
31 
(i 
201 
7 



a Extended source; see Sec. 3.5. 

k No SUMSS source within 30 arcsec, but SUMSS source and ATCA detection within 35 arcsec; see Sec. 6. 

c Associated with NGC 1824; see Sec. 3.5. 

ri Near deep SZ decrement but not spurious; see Sec. 6. 

c Offset between 1.4 mm and 2.0 mm emission; see Sec. 3.5. 

f Associated with NGC 1853; sec Sec. 3.5. 

s Possibly spurious detection from sidelobe of strong SZ decrement; see Sec. 6. 

Chance superposition (45 arcsec separation) of SUMSS source and new dust-dominated detection; see Sec. 3.5. 
1 Two SUMSS sources within 35 arcscc; see Sec. 6. 

j No SUMSS source within 1 arcmin, but associated with HE 0536-5435 [VCV2001]; see Sec. 6. 

k Associated with /5 Pictoris; see Sec. 6. 

1 Associated with ESQ 160- G 002; see Sec. 3.5. 



